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ABSTRACT 

Recent X-ray and weak-lensing observations of galaxy clusters have revealed that 
the hot gas does not always directly trace the dark matter within these systems. Such 
configurations are extremely interesting. They offer a new vista onto the complex in- 
terplay between gravity and baryonic physics, and may even be used as indicators 
of the clusters' dynamical state. In this paper, we undertake a study to determine 
what insight can be reliably gleaned from the comparison of the X-ray and the weak 
lensing mass maps of galaxy clusters. We do this by investigating the 2D substructure 
within three high-resolution cosmological simulations of galaxy clusters. Our main 
results focus on non-radiative gas dynamics, but we also consider the effects of ra- 
diative cooling at high redshift. For our analysis, we use a novel approach, based on 
unsharp-masking, to identify substructures in 2D surface mass density and X-ray sur- 
face brightness maps. At full resolution (~ l'bh~'^ kpc), this technique is capable of 
identifying almost all self-bound dark matter subhaloes with M > \Q'^^h~^Mcr). We 
also report a correlation between the mass of a subhalo and the area of its correspond- 
ing 2D detection; such a correlation, once calibrated, could provide a useful estimator 
for substructure mass. Comparing our 2D mass and X-ray substructures, we find a 
surprising number of cases where the matching fails: around one third of galaxy-sized 
substructures have no X-ray counterpart. Some interesting cases are also found at 
larger masses, in particular the cores of merging clusters where the situation can be 
complex. Finally, we degrade our mass maps to what is currently achievable with 
weak-lensing observations 100/i~"'^kpc at z — 0.2). While the completeness mass 
limit increases by around an order of magnitude, a mass-area correlation remains. Our 
paper clearly demonstrates that the next generation of lensing surveys should start to 
reveal a wealth of information on cluster substructure. 

Key words: X-rays; galaxies: clusters ~ methods: numerical - gravitational lensing. 



1 INTRODUCTION 

Tlie advent of tlie weak lensing technique has allowed ob- 
servers to directly probe the distribution of mass in galaxy 
clusters, rather than simply assuming that light provides 
an accurate tracer of the underlying dark matter (DM) dis- 
tribution. This allows us to separate out shortfalls in our 
understanding of baryonic physics from direct challenges to 
the Cold Dark Matter (CDM) model, providing an excellent 
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test of the predictions of the CDM paradigm itself and a 
clearer picture of the influence of the baryonic component. 

In recent years, there has been a flurry of papers, with 
these goals in mind, which compare weak lensing mass 
reconstructions to X-ray images of galaxy clusters (e 
ISmail et al.]|l997l : iMachacek et aLlbOOa ISmith et al1l200: 
Some such observations have highlighted dramatic excep- 
tions to the basic picture t hat light follows mass, most 
famously, the bullet cluster (|Clowe et al.ll2004l ) where the 
main peaks in the X-ray image are offset from those in 
the weak lensing mass reconstruction. There have been 
several follow-up theoretical studies of this unique sys- 
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tern (for example, Takizawa 



20061 : ISpringel fc Farraij boOTi : 



iMastropietro fc BurkerlJl200S ) which conclude that its main 



features can be reproduced well by idealized, non-radiative 
merger simulations suggesting the driving factor is ram- 
pressure. 

There have also been observations of clusters with fea- 
tures in the weak lensing map which are absent in the X - 
ray image. For ex a,mple, in MS1054-0321 (|jee ct al' 2005) 
and in Abell 1942 (|Erben et al.ll2000l : [Crav et al,.2001. '), for 
which several theories are put forward: chance alignments of 
background galaxies, galaxy clusters that have not yet virial- 
ized and so possess little X-ray emitting gas or substructure 
within the cluster that has somehow been stripped of its 
gas. Even more puzzling is the recent observation of Abell 
520 (|Mahdavi et al.ll2007l ). in which an X-ray peak with no 
corresponding mass concentration and a mass concentration 
with no galaxies are detected. This is postulated to be a re- 
sult of either a multiple body interaction, or the collision of 
weakly self-interacting DM during the merger event. Most 
recently is the obse rvation of another extreme merger event 
l|Bradac e t al."2008\ in which two clusters with M~ lO^M© 
are both displaced from the single peak in the X-ray emission 
suggesting even higher mass substructures could be seen to 
be 'dark'. 

On the galaxy-mass scale, studie s of X-ray observations 
of the hot haloes of elliptical galaxies /Machac ek et al.ll2006l ) 
exhibiting features characteristic of ram pressure stripping 
were carried out, suggesting we should expect to find galaxy- 
sized subhalo es that are dark in X-rays. However, a system- 
atic study bv lSun et all (|2007l ) found 60 per cent of galaxies 
brighter than 2Lt still retained small X-ray coronae, poten- 
tially indicating a more complex picture than just hydro- 
dynamics, involving the suppression of heat conduction and 
viscosity by magnetic fields. 

There have been many theoretical studies with the aim 
of understanding the global properties o f purely DM sub - 
structure. For example, the systematics (|Gao et aLll2004h . 



i pl( „ .^-^ 

evolution l|Gill et al.fr2004l : iReed et al.ll2005l), effe cts of the 



parent halo merger history ("Tavlor fc Babulll2004l ) and spa- 
tial distribution (Diomand et al. 2 0041 ') of subhalo popula- 
tions have all been studied in great depth. Attention is now 
also b eing paid to the fate of the gas in subhaloes. iHested 
l|2006l ) incorporated a hot halo component into an analyt- 
ical model of ram pressure stripping of galaxies in groups 
and clusters and found that most galaxies were readily 
stripped of the majority of this. Inspired by the first obser- 
vations of cold fronts in Chandra data (e.g. lMarkevitch et al.l 
L2OOO) some authors invoked separa tions between the hot gas 
and DM of either the main clu ster ([Ascasibar fc MarkevitchI 
I2OO6I ) or a merging subcluster (jTakizawa 20051 ) as a possible 
mechanism for their production. There were also many other 
complimentary studies into the fate of gas in subhaloes on 
the gr oup or cluster mass scale. For example, Bialck et al. 
l|2002t l report the ablation of gas away from the core of a 
merging subcluster's DM potential, in a cosmologica l simu- 
lation, resulting in adiabatic cooling and lHeinz et al.l (|2003l ) 
use idealized merger simul ations to study th is process in 
more detail. More recentlv jPo ole et al.l l|2006l l performed a 
suite of idealized cluster mergers and found gas in the both 
cores was often disrupted, leading to additional transient 
structures in the X-ray emitting gas. In order to specifically 
investigate the fate of hot gas in galaxies orbiting in groups 



and clusters iMcCarthv et al.l (|2008l ) studied a suite of hy- 
drodynamic simulations. They find the majority of the hot 
gas is stripped within a few gigayears but that around 30 
per cent is retained even after 10 gigayears. 

Much of this work on the gaseous component, uses 
simulations of idealized mergers in order to reproduce spe- 
cific observational features of galaxy cluster substructure. 
What is required now is a similar treatment to that afforded 
for DM subhaloes; a systematic study of the statistics of 
hot gas substructure in fully cosmological simulations. In- 
deed there have only been two studies of this kind already, 
(|Tormen et al.ll2004l : lD"olag et al.l |2008j): The former focusses 
on the time evolution of subhaloes in non-radiative simula- 
tions, while the latter examines how the overall distribution 
of subhalo masses and compositions differ, depending on the 
physics incorporated. There are two main issues still to ad- 
dress. Firstly, many of the interesting substructures seen in 
X-ray images of clusters (tidal tails, diffuse gas clouds etc) 
are omitted from simulation studies which simply identify 
substructure as hot gas bound to subhaloes. Secondly, ob- 
servationally we can only view the substructure in projec- 
tion; how does this relate to the substructure in 3D? Both 
of these issues can be addressed by undertaking an analysis 
of galaxy cluster substructure in 2D, allowing projection ef- 
fects to be quantified without restricting the analysis to the 
bound components. A comprehensive study in this area will 
help us to construct a framework within which to interpret 
the surprising results from comparisons between weak lens- 
ing and X-ray observations, of which there will undoubtedly 
be many more in the near future. 

In this paper, we use high resolution resimulations of 
three galaxy clusters to compare the substructure in the 
hot gas and DM components and examine what factors af- 
fect their similarity, or otherwise. We use a technique based 
on unsharp- masking to identify enhancements to the cluster 
background in maps of the X-ray surface brightness and to- 
tal surface mass density, providing us with catalogues of 2D 
substructures. Our aims are to understand the relationship 
between 3D DM subhaloes and our 2D total mass substruc- 
ture catalogues, including the contribution of 3D subhaloes 
that lie infront of or behind the cluster, yet within the map 
region. We wish to understand how these 2D mass sources 
then relate to substructures in the projected X-ray surface 
brightness, in order that we may place some constraints on 
the frequency of mismatches between substructure in the 
hot gas and DM and the mass scales at which these occur. 
Finally, we investigate how various selection and model pa- 
rameters influence these two relationships. 

The paper is structured as follows. In Section[2l the sim- 
ulation properties, selection of the cluster sample and gen- 
eration of the maps are outlined. The detection technique 
and properties of our 3D subhalo catalogues are included 
in Section |3] while Section |4] provides the same information 
for our 2D substructure catalogues. In Section [S] the results 
of a direct comparison between the 2D mass map substruc- 
tures and the 3D subhaloes are presented. We investigate 
the likelihood of finding a 2D X-ray counterpart for each 2D 
mass substructure in Section |6] and explore the effect that 
redshift, dynamical state, the inclusion of cooling and ob- 
servational noise have on this in Section [7] Section [7] also 
includes several case studies, to illustrate in more detail the 
fate of a 2D mass substructure's hot gas component when a 
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2D X-ray counterpart cannot be found. We provide a short 
summary of our results at the end of Sections 5, 6 and 7, 
should the reader wish to skip to the end of these sections. 
Finally, Section |8] outlines the main conclusions and impli- 
cations of this work. 



2 CLUSTER SIMULATIONS 

We use the resimulation technique to study the clusters with 
high resolution. Th ree clusters were selected from the larger 
sample studied by iGao et al] (|2004l ') and resimulated with 
gas u sing the publicly-available GADGET2 A^-body/SPH 
code (jSpringell l2005l ). A ACDM cosmological model was 
assumed, adopting the following values for key cosmolog- 
ical parameters: f2m = 0.3, Qa = 0.7, fib = 0.045, h = 
0.7, (78 = 0.9. The DM and gas particle masses in the high- 
resolution regions were set to mjark ~ 4.3 x 10* /i^^M© and 
"^gas = 7.7 X 10^ /i~^Mq respectively, within a comoving 
box-size of 479/i~^Mpc. The simulations were evolved from 
2 = 49 to 2 = 0, outputting 50 snapshots equally spaced in 
time. The Plummer gravitational softening length was fixed 
at e = 10/i~^kpc in the comoving frame until z = 1, after 
which its proper length (e = 5/i^^kpc) was fixed. 

For our main results, we have chosen not to incorporate 
the complicating effects of non-gravitational physics (partic- 
ularly radiative cooling and heating from galaxies), for two 
reasons. Firstly, we wish to investigate any differences be- 
tween the hot gas and DM in the simplest scenario, i.e. due 
to ram-pressure stripping and viscous heating of the gas. 
Secondly, a model that successfully reproduces the observed 
X-ray properties of galaxy clusters in detail does not yet ex- 
ist, and so only phenomenological heating models tend to be 
implemented in cluster simulations. Nevertheless, we include 
a limited analysis of the effects of non-gravitational physics 
on our results, namely allowing the gas to cool radiatively 
at high redshift, in Section [7] We defer a study of the addi- 
tional effects of heating from stars and active galactic nuclei 
to future work. 

2.1 Cluster identification and general properties 

To define the properties of each cluster within the simula- 
tion data, DM haloes were first identified using a Friends 
of Friends (FoF) algorithm, where the position of the most 
bound DM particle was taken as the centre. The Spherical 
Overdensity (SO) algorithm was then used to grow a sphere 
around this centre and determine rsoo, the radius contain- 
ing a total mean density, p — 500pcrit(2), where pcrit{z) is 
the critical density for a closed universe at redshift, z. This 
radius was chosen as it approximately corresponds to the 
upper limit of the extent of X-ray observations. The three 
clusters are labelled A, B and C respectively. 

At the end of each simulation {z = 0), the masses of 
clusters A, B and C within rsoo were M500 = [6, 3, 6] x 
lO^*/i~^M0, approximately corresponding to [1.4, 0.7, 1.4] x 
10® DM particles respectively. To determine the evolution 
of each cluster with redshift, candidate progenitors were se- 
lected by finding all FoF groups at the previous output, 
whose centres lie within O.Srsoo of the present cluster's cen- 
tre. We adopt a short (typically b = 0.05, but decreasing to 
b — 0.025 for problematic snapshots) FoF linking length to 



avoid the linking of two merging progenitors in close proxim- 
ity, which can lead to fluctuations in the centre from output 
to output. We then select the most massive object as the 
progenitor, except when there are several candidates with 
similar mass, in which case we choose the one that is clos- 
est (this only occurs for cluster C). Our choice meant that 
where cluster C undergoes an almost equal-mass merger at 
z = 0.5, we did not end up following the most massive ob- 
ject at higher redshift, but tests confirm that this choice has 
no bearing on our main results. 

Fig. [1] (top panels) shows how M500 grows with time 
for each of the three clusters. For convenience, we have used 
redshift as the time axis; our study is limited to outputs 
between redshifts and 1. By design, the mass histories 
vary significantly between the clusters: cluster A undergoes 
several major mergers (leading to abrupt jumps in mass) 
early on, then settles down at z ~ 0.4; cluster B accretes 
matter over the whole period; and cluster C undergoes two 
massive mergers (at z ~ 0.9 and z ~ 0.4 respectively) with 
relatively quiet stages in between. 

2.2 Map making 

For our main analysis we constructed maps of surface mass 
density (dominated by the DM) and bolometric X-ray sur- 
face brightness for each cluster at each redshift. The former 
quantity is formally related to the volume density (p) as 

E = | pdl, (1) 

while the latter is related to both the electron density and 
temperature of the intracluster plasma 

^- = MTTIj^ /'^^^(^^)'^^' (2) 

although note that features are primarily due to fluctuations 
in the density. For the analysis that follows, the explicit 
redshift dependence of the surface brightness is ignored and 
we further assume the ideal conditions of inflnite signal to 
noise (except for discreteness noise due to the finite number 
of particles employed). The conversion from gas to electron 
density is performed assuming a fully ionised, Z = plasma 
(so Uc ~ 0.9pgas/ni-H) a nd the cooling rate is comput ed using 
the tables calculated bv lSutherland fc D opita (1993) for Z = 
O.SZq, the typical metallicity of the intracluster medium 
(ICM). 

The e stimation procedure followed is similar to that em- 
ployed by lOnuora et al.l (|2003D . Briefly, a cuboid is deflned, 
centred on the cluster, with sides of proper length 2r5oo in 
the X and Y directions and Sr^oo in the Z direction (to 
capture material associated with the cluster along the line- 
of-sight). Particles that reside within this cuboid are then 
identified and projected in the Z direction on to a 2D ar- 
ray of 400 X 400 pixels. The pixel size was chosen to sample 
length scales of at least the Plummer softening length (at 
z — 0, rsoo ~ l/i~^Mpc), so that all real structures were 
capable of being resolved by the map. 

The gas particles are not point-like, but spherical 
clouds of effective radius, h, and shape defined by the 
(sphne) kernel used by the ga dget2 SPH method (see 
ISpringel. Yoshida fc White! l200lD . Thus, all gas particles 
were smoothed onto the array using the projected version 
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Figure 1. Top panels: mass histories, Afsoo versus redshift, for cluster A (left), B (middle) and C (right). Middle panels: normalised 
RMS centroid displacement, a dimensionless estimator of a cluster's dynamical state (see text for details). The horizontal line indicates 
a value of 10 per cent (which we define as the threshold for a major merger) and filled squares in the top panels show outputs where the 
estimator exceeds this value. Bottom panels: examples of X-ray surface brightness maps for a cluster with low (0.01; left), intermediate 
(0.1; middle) and high (0.2) values of the RMS centroid shift, respectively. Here, contours illustrating equally spaced values of logarithmic 
surface brightness (white) and surface mass density (black) are overlaid. 



of the kernel. To reduce the noise in the mass maps, domi- 
nated by DM particles, densities and smoothing lengths were 
adopted in a similar way (defining h such that each DM par- 
ticle enclosed an additional 31 neighbours) and smoothed 
using the same kernel as with the gas. 

The projected mass density at the centre of each pixel, 
Ro, is therefore 

1 

E(Ro) = -— ^ m,w(|R, -RoUO, (3) 

^ I— 1 

where Apix is the pixel area, the sum runs over all TV particles 
within the cuboid region, rrii is the mass of particle i and w 



is the projected SPH kernel, suitably normalised to conserve 
the quantity being smoothed. 

The (redshift zero) X-ray surface brightness is estimated 
using a similar equation 

Ex (Ro ) = :r-r m.ne.i A(r,)m(lRi - Rol,/iO, 

^ 1—1 

(4) 

where the sum runs over all hot (Ti > lO^K) gas particles 
and we have assumed equivalence between the hot gas and 
electron temperatures. 

The maps are re-centred for analysis, such that the new 
centre is set to the brightest pixel in the X-ray surface bright- 
ness map, as would generally be the case for observations. 
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The allowed alteration is restricted to ensure that the centre 
does not 'jump' to a bright substructure (this is unrealistic, 
but possible because of our simple non-radiative model) as 
this would undermine the effort directed at following the 
assembling structure. 

The bottom three rows of Fig. [2] illustrate surface mass 
density (left column) and surface brightness (right column) 
maps for cluster A, B and C ai z — 0. Qualitatively, the 
strongest features are clearly present in both maps, but there 
are some differences, notably the lack of some obvious sub- 
haloes in the X-ray maps and extended features in the X-ray 
maps due to stripped gas. It is also clear that the brightest 
X-ray substructures tend to be much rounder than in the 
mass maps, as expected, since the gas traces the potential, 
which is smoother and more spherical than the density. 



2.3 Characterising dynamical state from the maps 

Our first application of the maps is to estimate the dynam- 
ical state of the clu ster from its v i sual a ppearance. We em- 
ploy the method of lO'Hara et al.l (|2006l ). also f ound to pro- 
vide a reliable indicator of dynamical state by iPoole et al.l 
l|2006t ). using idealised simulations of cluster mergers. For 
this method, the displacement between the X-ray peak and 
centroid is calculated within circular apertures ranging from 
rsoo down to O.OSrsoo in radius, decreasing by 5 per cent 
each time, and then the RMS value of the displacement is 
computed, relative to rsoo. We found this technique to be 
most effective when using heavily smoothed maps to com- 
pute the centroid, so adopt the smoothing kernel used in 
our substructure detection algorithm (described in Section 
|4|, but with a = O.lrsoQ. The position of the peak is always 
taken as the centre of the cluster, as defined in the previous 
subsection. 

The variation of this RMS 2D statistic with redshift is 
shown in the middle panels of Fig. [T] Values above 10 per 
cent (indicated with a horizontal line) are typically found 
whe n a cluster is unde rgoing a major merger (see, for exam- 
ple, [Poolele£|al][200^. The redshifts at which this value is 
exceeded are also indicated with filled symbols in the mass 
histories (top panels). The bottom panels of the same fig- 
ure give examples of clusters with low, intermediate and high 
values of the RMS centroid shift, clearly showing an increase 
in dynamical activity. 



3 3D SUBHALO DETECTION 

Although the key objective of our analysis is to study the 
X-ray and mass maps of the clusters, we can draw impor- 
tant insight from an analysis of the full 3D data. In this 
section we identify 3D self-bound DM subhaloes in the map 
region (a cylinder of radius rsoo and length Srsoo, centred 
on the main cluster) and investigate their global properties. 
The information we glean from this analysis will help us to 
interpret our results in Section [S] by allowing us to distin- 
guish the underlying physical mechanisms from any effects 
introduced by our method. 



3.1 Detection technique 

We use a version of SUBFIND (|Springel et al.llioOll ') to decom- 
pose FoF groups (identified for this purpose with b = 0.2) 
into 3D self-bound subhaloes. The modi fied ve rsion, kindly 
provided by Volker Springel (see also D olag et al.. ;2008| ), 
identifies both gas and DM particles (and star particles when 
relevant) within each subhalo. A region larger than the final 
map region is analysed such that all subhaloes that con- 
tribute significantly to, but may not be fully within, the 
map region are included. 

We employ a threshold of 100 DM particles, correspond- 
ing to a mass, M ~ 4 x lO^°/i~^M0, as our minimum reso- 
lution limit for the subhalo catalogues. As we will show, this 
is significantly below our 2D completeness limit (determined 
in Section [Sjl . As our 3D subhalo catalogues will form the 
basis for comparison with 2D substructure, we consider not 
only subhaloes that lie entirely within the map region, but 
also those with at least 75 per cent of their mass along the 
line of sight (as defined in Section [21), even if their centre 
co-ordinates are outside rsoo in projection. Note that, even 
if less than 100 per cent of the subhalo's particles are within 
the map region, the whole DM mass of the subhalo is still 
recorded. 

The mass of each subhalo is computed using only the 
DM particles, to reduce its dependence on the amount of gas 
stripping that has occurred (the mass, M^ub, therefore refers 
to the DM mass of the subhalo). We take the centre of the 
subhalo to be the position of the most bound particle, but 
also calculate a projected centre, to facilitate matching with 
the substructures in the map. This was determined to be 
the position of the peak projected number of DM particles, 
i.e. the co-ordinates of the cell containing the most particles 
when each subhalo's DM particles are binned according to 
their X and Y co-ordinates (particles with Z co-ordinates 
outside the map region are excluded). 

3.2 Properties of 3D subhaloes 

Before we begin discussion of the results in this section, it is 
important to note that we always include the main cluster 
halo in the subhalo data. This is important to facilitate the 
comparison to 2D substructures detected in the maps later 
on, as the cores of the clusters (see the central mass density 
peaks clearly evident in the first column of Fig. [5J are de- 
tected in 2D and these cores, therefore, are detected as 2D 
substructures in their own right. 

In Fig. [3] we show the cumulative subhalo DM mass 
function for subhaloes with their most bound particle in- 
side 3D rsoo, down to our imposed resolution limit of 100 
DM particles. The results for cluster A (solid), B (dotted) 
and C (dashed) are shown individually for z = 1 (left), 0.5 
(middle) and (right) in the first row. Note that the main 
cluster itself is the most massive subhalo. The total number 
of resolved subhaloes, ranging from less than 10 to nearly 
60, depends on cluster mass. For example, cluster C has sig- 
nificantly fewer subhaloes at z = 1 and z — 0.5 than the 
other two clusters, but has more ai z — 0. This increase re- 
flects cluster C's major merger at z ~ 0.4, as seen in Fig. [1] 
However, when the subhalo masses are scaled to the parent 
cluster mass, the scatter between clusters and redshifts is 
much smaller, as shown in the second row. 




Figure 2. Examples of cluster maps. Left column: logarithmic surface mass density maps for, top to bottom, cluster A cooling run, 
cluster A, B and C non-radiative runs at 2 = 0. Right column: the same but for the logarithmic X-ray surface brightness maps. 
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Figure 3. Subhalo mass functions for cluster A (solid lines), B (dotted lines) and C (dashed lines), at redshifts, z = I (left), z = 0.5 
(middle) and z = (right). First row: Cumulative subhalo (DM) mass functions. Second row: As for first row but with the subhalo DM 
masses scaled to Afsoo. Data are for subhaloes with their most bound particle within r^og. 
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Figure 4. Left panel: fraction of resolved subhaloes within the full map region that lie inside rsoo (solid lines) and 2r5oo (dashed 
lines). Dot-dash line represents fraction within rsoo expected for a uniform distribution (from ratio of volumes). Middle panel: fraction 
of subhaloes within rsoo (solid) and the the full map region (dashed) that have no hot gas (/gas = 0). Note that all lines are at 
zero in high mass bin. Right panel: as in the middle panel but for subhaloes that have little hot gas (/gas ^ 0.5 Ot/fim). Note that 
both black lines are at zero in high mass bin. In all panels, results for the same subhalo DM mass (10^^ < Mg^^,//^~^MQ < 10^^; 
lO^^ < Mgub//i~^M0 < 10^^; Msub > lOl^?i~^M0) and redshift (red lines: < z ^ 0.2; black lines: 0.5 ^ 2 < 1.0) intervals are shown. 
Note that data for all three clusters are stacked here and that the main clusters themselves are included in these data. 



We have also examined how the properties of subhaloes 
in the map region vary depending on whether or not they 
lie within rsoo in 3D, to assess the impact of subhaloes pro- 
jected along the line of sight. It is particularly important 
that we examine the distribution of subhaloes, because of 
the unusual geometry we are using (a cylinder, rather than 
a sphere). In the left panel of Fig. |4l we show the fraction 
of subhaloes (including the main cluster halo) in the entire 



map region that lie within rsoo (sohd) and 2r5oo (dashed) 
in 3D for the redshift intervals ^ z ^ 0.2 (red lines) and 
0.5 ^ 2 < 1.0 (black lines). These redshift intervals were cho- 
sen to include an equal number of snapshots (11). Around 
half of the low mass subhaloes lie within rsoo which is signif- 
icantly higher than for a uniform distribution, for which we 
would expect, — |vrr5oo/87rr|oo ~ 1/6 (indicated 

with the dot-dash line). Nearly all subhaloes (80 — 100 per 
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cent) lie within 2r5oo , suggesting that the contribution to the 
map from substructure outside the cluster's virial radius is 
small. The rise, compared to lower mass bins, in the fraction 
of subhaloes with Mgub > lO^^/i'^M© within rsoo is due to 
the presence of the cluster cores. The cluster cores dominate 
this bin (in number) and since the map region is centred on 
them, they are always within rsoo by design. At high red- 
shift, the fraction of galaxy-sized (Afsub < IQ^^hr^y^Q) sub- 
haloes within rsoo is approximately 10 per cent higher than 
at low redshift. This is likely to be caused by the effects 
of tidal forces, stripping the DM as the subhalo orbits in 
the cluster potential. This effect may reduce the likelihood 
of finding subhaloes which are dark in X-rays in this mass 
range at low redshift, since they may move to lower DM 
mass bins (via tidal stripping) shortly after their hot gas is 
removed. 

Given the aims of this investigation, we want to try 
to place some limits on the fraction of DM substructures 
without X-ray emission we expect to find. In the middle 
and right panels of Fig. |4l we now plot the fraction of sub- 
haloes (within rsoo (solid) and the full map region (dashed)), 
with no hot (T> lif K) gas (/gas = 0) and little hot gas 
(./gas ^ 0.5f2b/i^m), respectively. This somewhat arbitrary 
threshold, /gas ^ O.Silb/f^m, was chosen simply to distin- 
guish hot gas-poor subhaloes from hot gas-rich subhaloes. 
Note that, by definition, the rest of the subhaloes (1 — /subs) 
fall into the latter category and have fgas > 0.5r2b/^^m- 
The main trend apparent is that the fraction of empty or 
low- gas subhaloes is hig her at lower mass, in agreement 
with lTormen et ahl l|2004l ). who find the survival time of hot 
gas in subhaloes is a strong increasing function of subhalo 
mass. Without the added effects of radiative cooling and 
energy injection from galactic winds, for example, our re- 
sults already predict that the vast majority of galaxy-sized 
(10" < A/sub//i~^M0 < 10^^) subhaloes are substantially 
depleted of hot gas, while the opposite is true on group (and 
cluster) scales. We also find that more subhaloes have no hot 
gas at low r e dshift than at high redshift, in agreement with 
iDolag et ai] l|2008l ). We note that the middle panel of Fig.g] 
is insensitive to the temperature threshold, since the vast 
majority of subhaloes with no hot (T > 10^ K) gas have no 
gas of any temperature. 

The vast difference between the fraction of empty sub- 
haloes in the lowest mass bin and at higher masses (e.g. ~ 80 
per cent of subhaloes with lO" ^ Msub/ft"^Mo < 10^^ al- 
ready gas- free at high redshift, yet still only ~ 30 per cent 
with 10^^ ^ Msub//i~^M0 < 10 " gas-free at low reds hift) is 
qualitatively in agreement with iTormen et al] (|2004 ). if we 
assume higher redshift to indicate less time since infall. They 
find complete removal of hot gas within 1 gigayear (typi- 
cally massive galaxies) to 3 gigayears (typically groups) of 
entering the cluster's virial radius on average. Results from 
[McCarthv et al. (2008) are in general agreement, but indi- 
cate ~ 30 per cent of hot gas in a halo (typically a mas- 
sive galaxy) can survive much longer (~ 10 gigayears); a 
result shown to i mprove colours o f satellite galaxies in semi- 
analytic models (|Font et al.l2008l ) . We find that the majority 
of subhaloes with M > lO^'^/i~^M0 always retain some hot 
gas and indeed at least 20 per cent have /gas > 0.5f2b/r2m. 
This shows our results are compatible with subhaloes retain- 
ing some of their original hot gas, although it seems in most 
cases the majority is removed. 



iDolag et al.l l|2008l ) find that stripping is very efficient 
with ~ 99 per cent of all subhaloes in rvir being gas-free at 
2 = 0. Note that this percentage will be dominated by their 
low mass subhaloes which are most numerous (and most gas- 
deficient) and so compares well with the percentage (90 per 
cent) that we find in our low mass bin. It remains to be seen 
how much gas has to be stripped before the likelihood of 
detecting the substructure in both the total mass and X-ray 
surface brightness maps is affected. 



4 2D SUBSTRUCTURE DETECTION 

A number of authors have used 2D weak lensing maps 
and X-ray images of clusters, both to compare the spatial 
distr ibution of hot gas and dark matter in these objects 
(e.g. [ciowe et ah I l2004l : iMahdavi et al. I l2007l : iBradac et al.l 
[2000) and to help infer their dynamical state, by measur- 
ing the offset betw een the centres of these two components 
(|Smith et al.ll200"5l ). The scope of the information about the 
underlying 3D system which such 2D comparisons could po- 
tentially provide, has not yet been explored and the present 
study is the first attempt to do this. 

The key features of this piece of work are a simple, yet 
effective, technique for identifying substructure in 2D maps 
of simulated clusters, in combination with an easy-to-use 
method for mapping 2D mass substructures to both 2D X- 
ray substructures and 3D subhaloes. First, we analyse our 
'perfect' observations (i.e. noise-free maps). This allows us 
to establish how many projected mass and X-ray substruc- 
tures can, in principle, be uniquely identified despite projec- 
tion effects and the intensity of the cluster background. This 
approach also provides insight into the fate of a projected 
mass substructure's hot gas when an X-ray counterpart is 
not found; the maps (unlike 3D data) allow immediate vi- 
sual follow up and reveal interesting features of the stripping 
process. We will explore how much of this is observable with 
current techniques in Section [T]?] by degrading the map res- 
olution and adding noise to both map types. 

4.1 Detection technique 

The first step towards detection is to enhance the substruc- 
ture in the maps. For this purpose we use a method based 
on the unsharp-masking technique, in order to remove the 
cluster background. The unsharp-masking technique itself 
has already been used as a visual aid by highlighting small- 
sca le structure in X-ray image s of galaxy clusters, for exam- 
ple iFabian et al] l|2003l. |2005| ). The main advantage is that 
it does not rely on the cluster being circularly symmetric, 
recovering the distribution of substructure well even in the 
most complex scenarios (i.e. multiple mergers, as is some- 
times the case in our simulations, especially at high red- 
shift). 

The first stage of the procedure is to smooth the maps 
with a preliminary Gaussian filter. This could be used to 
emulate the point-spread-function of a real instrument, but 
here we set the Full- Width-Half-Maximum (FWHM) to sim- 
ply match the spatial resolution of the simulation, as our 
results are presented in the limit of no added noise (other 
than intrinsic discreteness noise due to the finite number 
of particles employed). Our maps contain a fixed number 
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(a) Pre-smoothed mass map. (b) Unsharpmask for mass map. (c) Unsharpmasked mass map (d) Substructures detected in 

after a cut has been made. mass map. 




(e) Pre-smoothed X-ray map. (f ) Unsharpmask for X-ray map. 



(g) Unsharpmasked X-ray map 
after a cut has been made. 



(h) Substructures detected in X- 
ray map. 



Figure 5. An example of our substructure detection procedure for cluster C at z = 0.162. The top row corresponds to the surface mass 
density maps and the bottom row to the X-ray maps surface brightness maps. See text for further details. 



of pixels (200) across rsoo, corresponding to a length scale 
for each pixel of around 5/i~^kpc at 2: = 0, which is the 
equivalent Plummer softening length of our simulation (held 
fixed in proper units over the redshift range of interest 
here). The minimum length scale that should be trusted is 
around 3 times this, corresponding to the extent at which 
the gravitational force law becomes perfectly Newtonian in 
the GADGET2 code. Furthermore, rsoo is smaller at higher 
redshift, so our pixel scale is also smaller. We therefore set 
FWHM = 15/i"^kpc (physical) for all maps studied in this 
paper. It should also be noted that the maps are generated 
to be larger in X and Y than required so that the larger 
maps can be analysed to avoid edge effects in the region of 
interest. Panels (a) and (e) in Fig. [5] illustrate examples of 
these pre-smoothed maps. 

The second stage is to convolve the pre-smoothed maps 
again with a broader Gaussian kernel, to create the unsharp- 
mask image, shown in Fig. O panels (b) and (f). Here we 
fix (T2 to be 0.05 rsoo (corresponding to a FWHM rang- 
ing from approximately 35/i~^kpc to 120/i~^kpc over the 
redshift range), which was deemed to be the most effective 
value from extensive testing. This twice-smoothed version 
of the map is then subtracted from the pre-smoothed map, 
leaving a map showing just the enhancements to the cluster 
background. 

Utilising the commutative, distributive and associative 
properties of convolution, it is possible to derive one function 
that, when convolved with the map image, produces the 
same result as the series of operations described above. The 
kernel used to generate the pre-smoothed map approximates 



the Gaussian function, which is given by 



Gprelim 

where the normalisation, 



iVie 



and, in this case. 
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which is set by the spatial resolution of the simulation. Sim- 
ilarly, the combined operations of pre-smoothing and gener- 
ating the unsharp-mask image is simply 



Gu.s. =iVi,2e ^("i+'P, 
where the normalisation is now 

iVl,2 = 



1 



(8) 



(9) 



and (T2 is 0.05rsoo- The function representing the entire pro- 
cedure. 



F = G. 
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Gr 



(10) 



which is a close approximation t o the Mexican -hat function 
or the matched filter defined bv iBabul (|l990l ). is shown in 
Fig.l. 

The size of substructures that are detected are depen- 
dent on the combination of the standard deviations of the 
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Figure 6. ID visualisation of the convolution kernel equivalent to 
the whole background-subtraction procedure. For this example, a 
typical value of 1 Mpc is used for the cluster radius, giving 
(72 = 50 kpc, (Ti ~ 6 kpc (as defined in equation [T]! and 
ro ~ 12 kpc (as defined in the ID analogue of equation IllH . 
Note that this is for illustration only and as such the area under 
each Gaussian is normalised to 1 between ±oo in ID before the 
difference of the two is taken. 



Gaussians used to obtain the final image. We derive an ex- 
pression that characterises the width of the kernel and there- 
fore the scale of substructure to which our technique is sen- 
sitive. The characteristic width of the function in Fig. [6] can 
be determined by calculating the radius at which the ampli- 
tude of the function is zero. The radius of the zero-points, 
ro, is given by, 



2 I „2 



(11) 



Since a2 is expressed in units of rsoo , it has a slight redshift 
dependency (e.g. for cluster B, rsoo ~ 0.34 /i~^Mpc at z = 1 
and 0.78 h'^Mpc at z = 0), meaning that more extended 
substructures will be detected at lower redshifts. However, 
the increase in the value of the kernel width is only of the 
order of 20 per cent of its maximum value over the range 
of redshifts studied, ^ z < I (e.g. ro = 0.014 h'^Mpc for 
z = 1 and 0.018 /i~^Mpc for z — 0, averaged over the 3 clus- 
ters). We are limited to detect only 2D mass substructures 
of the order of the size of the kernel and these 2D mass sub- 
structures will, of course, be associated with a 3D subhalo 
mass. Since the typical size of a 3D subhalo above a given 
mass is larger at lower redshift, due to the decrease in the 
critical density as the universe expands, the trend for the 
kernel to also be larger at lower redshift actually reduces 
the redshift-dependence of the minimum 3D subhalo mass 
which we can detect in 2D. 

In order to pick out the true substructures from other 
fluctuations, any pixels with values less than /i -I- Xa (where 
X is an integer, representing our detection threshold, and ^ 
and a are the mean and standard deviations of the residual 
substructure maps) are discarded. Examples of the resulting 
maps at this stage of the procedure are shown in Fig. [5] pan- 



els (c) and (g). Substructures are then deflned similarly to 
the FoF technique but in 2D, grouping together neighbour- 
ing pixels with values greater than the background level. 
Ellipses are fitted to these pixel groups by finding the eigen- 
vectors (corresponding to the direction of the semi-major 
and semi-minor axes) and the eigenvalues (whose square 
roots correspond to the magnitude of the semi-axes) of the 
moment of inertia tensor. This allows us to determine the 
extent, orientation and circularity of the 2D substructures 
(see below). 

We investigate three values of X for the projected total 
mass map: 1, 3 and 5, and evaluate which is most successful 
when the comparison with the 3D subhalo data is made in 
Section [5] It was found that the X-ray surface brightness 
maps respond slightly better to our technique due to the 
fact the gas distribution is far smoother (because it traces 
the gravitational potential) and contains fewer small-scale 
fluctuations, meaning that less stringent cuts are required in 
order to achieve the same results (upon visual inspection). 
Therefore, the selection of the parameters used to define the 
catalogue of X-ray substructures is undertaken separately 
to that for the mass substructures. We found that X — 5 
is too stringent for the X-ray maps, removing substructures 
that are clearly visible by eye, whereas X = 1 and X = 3 
produce reasonable results for both the X-ray and the mass 
maps. 

4.2 Properties of 2D substructures 

The total number of substructures detected in the 1, 3 and 
5(T total mass catalogues (which consist of 90 maps, i.e. 30 
per cluster) is 3224, 1233 and 680 respectively. It is clear 
from these numbers that, as would be expected, the higher 
the value of X, the lower the number of detections. There 
are also considerably fewer X-ray substructures than total 
mass substructures for the same X value, with 1169 in the 
la X-ray catalogue and 707 in the 3a" X-ray catalogue. This 
can in part be attributed to the smoother distribution of the 
hot gas, which responds differently to the unsharp-masking 
procedure. However, it is also apparent (on visual inspec- 
tion) that there is simply less inherent substructure in the 
X-ray maps, particularly on small scales. 

First, we examine the distribution of total mass (solid) 
and X-ray surface brightness (dashed) substructure areas, 
^Bub, in the left-hand panel of Fig. [T] Aaub is defined as the 
number of pixels attributed to the 2D substructure in the 
unsharp- masked image multiplied by the physical area of the 
individual pixels. The distributions are very similar, except 
the X-ray surface brightness distribution peaks at a slightly 
higher value of Aaub , suggesting the X-ray substructures are 
typically more extended (this is confirmed by visual inspec- 
tion). There is little dependence on choice of catalogue here. 

We examine the shape of the substructures by plotting 
the distribution of circularities in the middle panel of Fig. [7] 
(line styles as before). Here, we define circularity, c = (b/a) 
where a is the major axis and b the minor axis of the ellipse. 
This distribution is very stable to choice of catalogue, sug- 
gesting that the morphology of the objects we detect changes 
little between catalogues. The distribution of mass substruc- 
tures peaks at c ~ 0.75, due to the triaxial nature of the 
DM su bstructure; the gas is slightly rounder and peaks at 
c~ 0.8.|K^ibi et al.l (|2008l ) obtain a similar result for their 
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Figure 7. Left: Area distribution of substructures. Middle: Distribution of circularities for substructures. Right: Angle between semi- 
major axis vector of substructure and radial vector from substructure centre to cluster centre. Surface mass density (solid) and X-ray 
surface brightness (dashed) 3cr catalogues. Vertical lines indicate the mean values in each case. 



projected sphericity of DM subhaloes, computed from the 
particles directly. This indicates that our detection method 
recovers the true 2D shape of the substructures successfully. 

The right-hand panel of Fig. [7] shows the distribution 
of the radial alignment of the 2D substructures with respect 
to the cluster centre. This is computed by first calculating 
the angle of inclination between the cluster centre (in 2D) 
and the centre of the ellipse representing a substructure. 
The alignment, 8, is then found by subtracting from this 
the angle of inclination of the semi-major axis of the ellipse. 
The range of 8 can be reduced to — 90° by treating oppo- 
site directions of the semi-major axis vector as equivalent. 
A mild tendency towards alignment is exhibited by the total 
mass substructures (solid), however the X-ray substructures 
(dashed) show no preferred direction. Knebe ot al. (2008) 
perform a similar calculation for their projected DM sub- 
haloes and found a much stronger tendency for alignment 
than we see here, when they considered all particles associ- 
ated with the subhalo. However, they investigate the effect 
of varying the percentage of particles they analyse by limit- 
ing the alignment measurement to the inner regions of the 
subhaloes. The trend for alignment shown in their results 
weakens as smaller percentages of particles are considered 
and comes into agreement with observations when analysing 
the inner 10 — 20 per cent. Our result is also in much better 
agreement with theirs for this region. This reflects the fact 
that our 2D detection technique finds only the cores of the 
original 3D subhaloes, which is demonstrated by the small 
scale of the detected 2D substructures. Therefore, we are ef- 
fectively performing our alignment and circularity analysis 
on o nly the innermost p articles and so find best agreement 
with iKnebe et al.l (|2008l ) when they similarly restrict their 
analysis. 



5 COMPARISON OF 2D MASS 

SUBSTRUCTURES TO 3D SUBHALOES 

In this section, by comparing the 2D total mass substruc- 
tures (described in Section |4]) with 3D self-bound DM sub- 
haloes (described in Section |3]), we assess the reliability of 
our 2D detection method and infer the 3D properties (e.g. 
subhalo mass) of our 2D substructures. We examine the 
completeness (with respect to 3D) of our 2D catalogues. 



as well as the number of individually resolved high-mass 
objects they contain, in order to select one total mass cat- 
alogue that is most suited for the analysis in later sections. 
Our catalogues contain substructures identified in all three 
clusters and at all redshifts (0 ^5 2: < 1). 

Ideally, we want to be 100 per cent complete down to at 
least Msub ~ lO^^/i~^M0, as this is the typical mass scale 
of substructures detected in current observations of the un- 
usual systems discussed in Section [T] However, high com- 
pleteness at lower masses would be desirable as smaller sub- 
halo es are the more likel y ones to be found stripped of their 
gas l|Tormen et al.l [20041 ). An additional criteria we wish to 
place on any detections is that, ideally, they are individu- 
ally resolved (i.e. not confused with another subhalo that 
is nearby in projection). We also look at the purity of our 
2D substructure catalogues by assessing the fraction of 2D 
substructures which we fail to associate with 3D subhaloes. 



5.1 Completeness 

Firstly, we determine the completeness of each of our three 
2D mass substructure catalogues {X = 1, 3, 5). This is done 
by starting with the 3D subhaloes and looking for 2D coun- 
terparts in the mass maps, then examining the resulting 
matching success (i.e. the fraction of 3D subhaloes for which 
a 2D counterpart is found) per 3D subhalo mass bin. The 
criterion for matching the 3D subhaloes to the 2D mass sub- 
structures is that the centre of the 3D subhalo must lie 
within the ellipse that characterises the 2D substructure 
(with a ~ 20 per cent margin for error, which was deter- 
mined by experiment). As discussed in Section|31 the default 
3D centre is taken to be the (projected) position of the most- 
bound particle in the subhalo, as identified by SUBFIND. This 
is a robust choice, comparing very well with the peak surface 
density in the maps in the vast majority of cases. However, 
during a complicated merger, we found that the most-bound 
particle can occasionally lie outside the cluster core (see Sec- 
tion 7.2), in which case we apply the position of the peak 
projected DM particle density of the cluster instead. 

Multiple 3D subhaloes can be matched to the same 
substructure in the mass map; we refer to this as a mul- 
tiple match. However, 2D substructures cannot share a 3D 
subhalo as our criterion means each subhalo is only ever 
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Figure 8. Fractional matching success per mass bin of 3D sub- 
haloes to 2D mass substructures as a function of subhalo DM 
mass for Icr 2D catalogue (solid line), 3(t 2D catalogue (dotted 
line) and 5(7 2D catalogue (dashed line). Bins are equally spaced 
in log(Msub)- 



Figure 9. Fraction of detected subhaloes per mass fim which are 
obscured (see text for definitions) as a function of subhalo DM 
mass for la 2D catalogue (solid line), Scr 2D catalogue (dotted 
line) and 5a" 2D catalogue (dashed line). Bins are equally spaced 
in log(Afsub). 



matched to one 2D substructure. It should also be noted 
that we start with a limited 3D catalogue containing only 
those subhaloes whose centres are within the projected rsoo 
and match to the complete catalogue of 2D substructures, 
which extends slightly beyond the projected rsoo (i.e. out- 
side the map). This simply prevents the failure to match 
a genuine 2D-3D pair when one substructure's centre lies 
slightly outside this boundary. 

Fig. [5] illustrates the completeness of our 2D catalogues 
as a function of subhalo mass (note that the main clusters 
are included in these data). Specifically, it shows the frac- 
tion of all subhaloes in the map region that are detected, 
including those subhaloes associated with the same 2D sub- 
structure, due to source confusion or genuine projection ef- 
fects (detailed below). In all three catalogues we clearly as- 
sociate 2D substructures with all 3D subhaloes that have 
Msub > lO"/i"^M0. The la, 3cr and 5(7 catalo gues are 90 
per cent complete per mass bin down to 3 x 10^^, 10^^ and 
3 X lO^^/i"^M0 respectively. 

The cut-off in completeness, below which our ability 
to retrieve 3D subhaloes from the projected data decreases 
sharply with mass, is a result of several limiting factors: the 
map resolution (effectively set by the pre-smoothing kernel 
size, (7i), the choice of (72 and simply the intensity of the 
cluster background. Low mass subhaloes have poor contrast 
against the background since they add little mass in addition 
to the total mass along the line of sight and so are difficult 
to distinguish. The mass at which this cut-off occurs is most 
sensitive to a\. As we demonstrate in Section [7.41 when we 
increase a\ by a factor of around 10 (more typical of the 
resolution of weak lensing mass reconstructions), the 90 per 
cent (per mass bin) completeness limit for the 3a catalogue 
becomes ~ IO^^/i'^Mq (see Fig. [26)) . 



5.2 Projection and Confusion 

Visual inspection of the projected mass maps reveals that 
two peaks that are very nearby can be detected as one 2D 
substructure (i.e. confused), if the lower density pixels be- 
tween them are not removed when pixels < /i -f Xg are 
discarded. This is not a concern if the mass ratio of the 3D 
subhaloes that have given rise to the 2D peaks is high (as 
the inclusion of the less massive object has little effect), or 
if they are both low mass subhaloes (~ 10"/i"^Mo), be- 
low the mass range we are interested in. However if, for 
example, a subhalo with a mass, M,ub ~ 10"/i-^Mo and 
the main cluster core give rise to two adjacent peaks in the 
map which are confused as one 2D substructure, we limit 
our opportunities to study the properties of the subhalo in 
detail. This is particularly important since subhaloes with 
Msub > IQ^^hr^yiQ are relatively rare. A related effect is 
that of projection, where two subhaloes that are aligned 
along the line of sight give rise to only one peak in the pro- 
jected mass map. 

Here we do not distinguish explicitly between projec- 
tion and confusion. Instead, we define a detected subhalo as 
obscured if it is part of a multiple match and is not the most 
massive subhalo involved; we would not consider such a sub- 
halo to be individually resolved. Fig.[9]shows the fraction of 
detected subhaloes per mass bin which are obscured. As ex- 
pected the obscured fraction at the high-mass end is lowest 
for the 5a catalogue and highest for the la catalogue, since 
the former is the most stringent when removing low density 
pixels between adjacent substructures, allowing them to be 
individually resolved. The trend reverses at low mass, how- 
ever, because the removal of low density pixels also erases 
small 2D substructures. Since this is more effective with a 
larger value of sigma, the fraction of obscured substructures 
(detected only because of their association with larger sub- 
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Figure 10. Purity of our sample, i.e. tlie fraction per area bin 
of 2D mass substructures in our 3cr catalogue matched to 3D 
subhaloes as a function of the physical area of the substructure, 
j4-aub. The vertical line marks A^^h = 3 X 10~*/i~^Mpc^, above 
which our sample can be taken to be pure. Bins are equally spaced 
in log{Asub). 



Figure 11. Correlation between physical area, Agubi of 2D total 
mass substructure (Scr catalogue) in unsharp-masked image and 
DM mass of 3D subhalo, Afgub, to which it is matched. Filled 
squares show SUBFIND background haloes (see text for details). 
Vertical line shows purity threshold for 2D mass catalogue and 
the other is best-fitting line given by equation 1121 



structures) increases. For the Scr catalogue, around 70 per 
cent at 10"/i"^Mq, 5 per cent in the 10^^ - IO^^/i'^Mq 
mass range and zero at the high-mass end, are obscured. 
On inspection of the maps, it is apparent that the obscured 
fraction at mass scales of ~ lO^^/i~^M0 typically occurs in 
the final stages of a merger and results from confusion when 
the two objects coalesce. 

We adopt the 3a catalogue from now on as it offers a 
small reduction in the obscured fraction at high masses while 
maintaining good completeness above Msub = 10^^ h~^MQ, 
detecting 98 per cent of all subhaloes above this mass. 

5.3 Purity 

We now consider the purity of our Scr catalogue, by under- 
taking the matching procedure in reverse, i.e. starting with 
the 2D mass substructures and trying to identify SD subhalo 
counterparts for these. The matching success (i.e. the frac- 
tion of 2D mass substructures that are successfully matched 
to a 3D subhalo, in this case) then provides a measure of 
the purity. This is important as it tells us in which regions 
of parameter space the raw 2D data could potentially be 
used directly, without reference to the SD data for calibra- 
tion. Fig. [To] shows the fractional matching success of 2D 
mass substructures to SD subhaloes versus the characteris- 
tic physical area of the 2D substructure, As^b- 

We achieve very high purity down to A^ub — S x 
10~*/i~^Mpc^, close to the approximate projected area of 
our combined kernel, F, {nrQ ~ 10~"^/i~^Mpc^; see equa- 
tions [10] and [11] for definitions of F and a;o). Reasons for not 
finding a SD subhalo to match every 2D substructure are 
three-fold. Firstly, we have detected a substructure associ- 
ated with a SD subhalo with less than 100 DM particles (i.e. 
our minimum allowed subhalo size) . Secondly, the substruc- 



ture detection is 'false' i.e. we have detected a transient en- 
hancement which does not constitute a self-bound subhalo. 
Or finally, there is an associated subhalo but matching has 
failed (matching becomes increasingly difficult as substruc- 
tures become smaller). 



5.4 Correlation between mass and area 

In Fig. 111! we take all 2D total mass substructures in the 
Scr catalogue which have been matched to SD subhaloes and 
examine the correlation between the area of the 2D substruc- 
ture, Asub and the DM mass of its SD subhalo counterpart, 
Msub. The figure contains data for all three clusters at all 
redshifts in the range ^ z < 1. The filled squares indicate 
objects which are SUBFIND background haloes as opposed to 
subhaloes. Background haloes consist of the most massive 
subhalo found in each FoF group plus any additional group 
particles that are gravitationally bound to it and which do 
not already belong to a subhalo. Effectively, the background 
halo is the parent halo in which the other subhaloes reside. 
The background haloes grouped in the top right are the clus- 
ter cores themselves, forming a separate population because 
the 2D detection corresponds to the core only, whereas the 
mass is that of the entire cluster. The background haloes at 
lower masses are smaller parent haloes which lie in front of 
or behind the main clusters. This figure shows that when we 
are above the completeness limit in terms of associated SD 
subhalo DM mass (Msub = 10^^/i"^Mq), most 2D substruc- 
tures are also in the region where we know our 2D catalogue 
is pure (i.e. Asub > S x 10~*/i~^Mpc^); the converse does 
not hold, however. 

For the rest of the substructures, that are not back- 
ground haloes, we demonstrate a strong correlation between 
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the 3D DM subhalo mass (Msub) and the 2D area (Asub). A 
least squares fit to this correlation yields, 

(12) 

where all points with Asvih > 3 x 10 "^h ^Mpc^ were 
considered. Using this correlation we can select a new thresh- 
old of Asub = 10~^ft~^Mpc^ (corresponding to a mass of 
lO^^/i^^M©) producing a sample of substructures with both 
high purity and high completeness (we refer to such cata- 
logues as pure). 

We can also estimate the intrinsic scatter in this relation 
using, 



where Mi is the mass value of each data point and 
log (Mfit) is the value computed using equation 1121 for the 
corresponding area. We find f-iog (jv4„t,) ~ 0.35 which sug- 
gests that the typical uncertainty in the DM mass of a sub- 
halo is around a factor of 2. For comparison, the fit was also 
made using the discarded la and 5cr 2D mass catalogues in- 
stead and the intrinsic scatter in the resulting relations was 
very similar (0.30 and 0.37 respectively) suggesting the qual- 
ity of the fit is independent of catalogue choice. Furthermore, 
for given value of A^Mh, the maximum variation in the esti- 
mated value of Afsub when comparing all three 2D catalogues 
with each other is approximately a factor of 3, comparable 
to the error from intrinsic scatter. We also note that the in- 
trinsic scatter is greater at the higher redshift, for example 
it is 0.31 when fitting only to data for ^ z ^ 0.2 and 0.41 
for 0.5 ^ z < 1.0. Such a correlation, though calibration- 
dependent, is potentially useful for providing a quick, rough 
estimate of subhalo DM mass determined from the observed 
area of a substructure in a weak-lensing map (assuming the 
substructure is resolved). 

5.5 Summary 

We have matched our 2D total mass substructure catalogues 
to self-bound 3D subhaloes and have identified the 3cr cata- 
logue as the most suitable for the analysis in future sections. 
This catalogue is at least 90 per cent complete in all subhalo 
mass bins above lO^'^/i~^M0 and pure above a projected 
area of 3 x 10~*/i~^Mpc^, which is close to the resolution 
limit of our kernel. We also note a strong correlation be- 
tween the (observable) area of the 2D substructure and the 
DM mass of the 3D subhalo. Using this we derive an area 
threshold, 10~^/i~^Mpc^, above which our substructure cat- 
alogues have both high purity and completeness. Projection 
and confusion effects above the completeness limit are min- 
imal. 



6 COMPARISON OF SUBSTRUCTURE IN 
THE HOT GAS AND DARK MATTER 
COMPONENTS 

We now address the main aim of the paper; comparing the 
substructure in the X-ray surface brightness and total mass 



maps. Again, we apply a simple matching technique, this 
time to our pairs of maps and then attempt to explore the 
underlying physical mechanisms which govern the resulting 
matching success, whilst also trying to constrain any poten- 
tial biases our method may have introduced. 

The catalogues of 2D mass substructures and 2D X-ray 
substructures are compared for each snapshot. The criterion 
for a match is that there is some overlap of the region en- 
closed by the ellipse that characterises the mass substructure 
and the region enclosed by the ellipse that characterises the 
X-ray substructure. In order to keep our method simple, we 
allow both 2D total mass and 2D X-ray surface brightness 
substructures to be matched to more than one substructure 
of the other type, rather than using additional matching cri- 
teria to prevent this. We use the term single match to refer 
to a unique pairing of one mass substructure with one X-ray 
substructure and the term multiple match for a mass sub- 
structure which has been matched to more than one X-ray 
substructure or vice versa. 

6.1 Direct matching 

An important feature of this work is the use of 2D data 
(maps) so, with this approach in mind, we first undertake 
the matching with no reference to the 3D subhalo data. This 
will allow us to confirm how reliable a picture the 2D data 
alone can provide as we can later compare our results to 
those which have been calibrated against the 3D subhalo 
information. 

As in Section[5]we undertake the matching procedure in 
two ways; starting with the 2D total mass substructures and 
seeking an X-ray counterpart for each and then starting with 
the 2D X-ray substructures and seeking a mass counterpart 
for each. Table [1] summarises the results of these matching 
processes, where the data for 2D mass substructures comes 
from the former and that for 2D X-ray substructures from 
the latter. Here we use the subscripts TM and SB to signify 
substructures in the total mass and X-ray surface brightness 
maps, respectively. 

Firstly, it is encouraging that the fraction of sub- 
structures which are matched to more than one object 
( /multimatchcd ) is very low (2 — 10 per cent) regardless of 
choice of 2D X-ray catalogue or whether only the pure sam- 
ple of 2D mass substructures is used. High numbers of single 
matches are preferred as this suggests the effect of confusion 
is limited and that the number of false matches is low. 

The ratios of X-ray to total mass substructures (■^^) 
show that when using the full 2D mass catalogue, there is 
always a dearth of X-ray substructures and so, regardless of 
the criterion employed, there will be unmatched 2D mass 
substructures. However, when using the pure 2D mass cata- 
logue, there is a factor of ~ 2 — 3 more X-ray substructures. 
The fraction of total mass substructures that are matched 
(/matched) roughly doublcs when moving to the pure sam- 
ple, suggesting that the majority of mass substructures dis- 
carded to obtain purity did not have an X-ray counterpart. 
This could be interpreted in one of two ways; the discarded 
2D mass 'substructures' were false detections and so one 
would not expect to find any corresponding substructure in 
the X-ray emitting gas, or they were real but may have cor- 
responded to low mass 3D subhaloes which are less likely to 
have retained their hot gas. In fact, around 75 per cent of 2D 
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Table 1. General results of direct matching between the 2D 
total mass and 2D X-ray surface brightness substructure cata- 
logues. Top rows: matching between all substructures in the X-ray 
catalogue, to all those in the total mass catalogue. Bottom rows: 
direct matching of all substructures in the X-ray catalogue to a 
pure xM > 10~^/i~^Mpc^) total mass catalogue. Columns: 

Multiple of (T used in indicated catalogue {X), Ratio of number of 
substructures ( -S^^ ) , fraction of total substructures in the cata- 
logue indicated that are successfully matched (/matched): fraction 
of all matched substructures in the catalogue indicated that are 
matched more than once (/multimatched)- 



Xtm ^se 



A^TM 



TM SB 



Imu It i matched 



TM SB 



Full 2D total mass substructure catalogue 
3 1 0.95 0.43 0.45 0.06 

3 3 0.57 0.32 0.59 0.08 



0.07 
0.03 



Pure 2D total mass substructure catalogue 

3 1 2.83 0.77 0.30 0.10 0.04 

3 3 1.71 0.67 0.43 0.10 0.01 



mass substructures below the purity threshold were matched 
to a 3D subhalo, suggesting it is the latter effect that domi- 
nates. Interestingly, when moving from the la to the 3a X- 
ray catalogue the fraction of mass substructures matched de- 
creases, but the same quantity for the X-ray increases. Here, 
the 3a X-ray catalogue is more pure as a greater fraction of 
its substructures can be matched to 2D mass substructures, 
however the la X-ray catalogue is more complete since a 
greater absolute number of its substructures are matched 
to 2D mass. A similar trade-off between purity and com- 
pleteness was seen when matching 2D mass substructures 
to 3D mass subhaloes in Section [5] The added complication 
here, of course, is that unlike the 2D mass substructures 
and 3D subhaloes, we cannot assume a 1:1 correspondence 
between the 2D mass and 2D X-ray substructures (in fact 
the deviation from this is the motivation for this work), so 
a completeness limit cannot really be established. 

Even more surprising than the large fraction of 2D mass 
substructures with no X-ray counterpart, is that a 2D total 
mass counterpart cannot be found for a high percentage of 
the 2D X-ray substructures. Even when considering the 3a 
X-ray catalogue, which picks out only the most defined 2D 
substructures in the hot gas, and matching this with the full 
3a 2D mass catalogue, 40 per cent of the X-ray substruc- 
tures still go unmatched. Investigating the properties of the 
matched and unmatched substructures should provide in- 
sight into this result. 



Focussing first on those substructures that are matched, 
we compare the properties of the 2D matched pairs. Fig. 
1121 demonstrates the tight correlation between the area of 
singly-matched X-ray (la catalogue) and total mass (3a 
catalogue) substructures. The best-fitting line for matched 
pairs from this combination of catalogues, where the 2D 
mass map substructure is in the pure region (-Asub,TM ^ 
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Figure 12. Correlation between areas of matched pairs of 2D to- 
tal mass, Agub tMi a^nd 2D X-ray, A^^^, sbi substructures. Data 
for single matches only (see text for definition) for matching be- 
tween 3f7 total mass and lo X-ray catalogues (data set using 3a 
X-ray catalogue has been omitted for clarity but shows a sim- 
ilar distribution). Solid lines are a least squares fit to all the 
data (extended line) and just that above the purity threshold 
(j^sub TM ~ 10~^/i~^Mpc'^; indicated with a vertical line) for 
the 2D mass catalogue; The latter is given by equation 1141 The 
complete data set has a correlation coefficient of 0.78. 



log ,^1^^^"^ = (0.83 ± 0.04) log ^;"'''™^ - (0.4 ± 0.1) 



/!.-2Mpc^ 



/i-2Mpc^ 



lO""^/!" Mpc"^) is given by. 



(14) 

Including all multiple matches as well increases the scat- 
ter, but the relationship is still clearly evident. Note that the 
X-ray substructures are slightly larger than the total mass 
substructures; this is partly due to the use of the la cat- 
alogue here, but also due to the more extended nature of 
the hot gas (for comparison, the gradient when using the 
3a X-ray catalogue is 0.91 ± 0.06, still less than 1). The 
outliers above the line can mostly be attributed to small 
2D mass substructures being matched to highly elliptical 
2D X-ray substructures, which are usually features near the 
centre of the main cluster corresponding to subhaloes ac- 
tively undergoing stripping. In many cases it is impossible 
to tell whether the match is valid or not, however the scatter 
occurs below the threshold area which demarks where our 
catalogue is pure (Aaub,TM ~ 10~^/i~^Mpc^). With this in 
mind it is unsurprising that some of the scatter here also re- 
sults from spurious detections, i.e. mass substructures that 
are later found not to correspond to a subhalo. Scatter be- 
low the line seems to arise from two situations 1) a small 
gas feature is detected that overlaps with a large mass sub- 
structure which has been stripped of its gas, i.e. the two are 
in chance alignment, 2) the match appears genuine yet the 
gas substructure is small, suggesting the outer regions of gas 
have already been stripped. 

Fig. [T3] shows the fraction of all 2D total mass sub- 
structures matched to X-ray (top) and the fraction of all 
2D X-ray substructures matched to total mass (bottom) as 
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a function of substructure area. Above 10"'^ Mpc^ the 
matching success is > 50 per cent per bin for both cata- 
logue types, however this still suggests a very high num- 
ber of substructures do not have counterparts in the other 
map. There are also many large unmatched substructures, 
for example around 10 per cent of 2D mass substructures 



in the -2.25 < log (Asub/ft" Mpc^ 



-2.0 range. Using 



equation [12] we can infer that this corresponds to a mass 
of approximately IO^'^/i'^Mq, suggesting these correspond 
to fairly massive 3D subhaloes. Increasing the radius of the 
unsharp-masking kernel used to detect the X-ray substruc- 
tures, to twice that of the fiducial kernel (i.e. O.lrsoo), yields 
a similar matching success. This indicates that the results 
of the substructure comparison that are shown here, do not 
depend significantly on this aspect of the substructure de- 
tection procedure. 

It is interesting that there is also a significant number 
of unmatched X-ray substructures, even at large areas. One 
would initially assume that once hot gas is separated from 
its DM subhalo it would disperse and so not be detected as 
a stand-alone substructure. Large unmatched 2D X-ray sub- 
structures were followed up individually by visual inspection 
of the maps and it appears that there are three main cate- 
gories. These are: 1) clearly defined substructures which are 
so displaced that they cannot visually be associated with one 
particular DM substructure (although there are typically 
candidates in the vicinity), 2) clearly defined substructures 
that are slightly offset from a nearby dark mater substruc- 
ture, and 3) detections of gas 'features' in the vicinity of 
the core during merger events - these 'substructures' cannot 
be directly associated with a DM substructure and it is not 
necessarily appropriate to do so. 

Scenario 1 incorporates the most clear-cut examples of 
2D X-ray substructures which are indisputably unmatched, 
whereas scenario 2 also includes those whose definition as 
matched or unmatched is somewhat subjective, as it is clear 
which mass substructure they belong to even though they 
are spatially distinct from it (for our purposes we call these 
unmatched). An example of the displacement of the X-ray 
component of a substructure can be seen in Section 17.21 
(see Case Study 2, Fig. I19|> . This example is a simple one, 
however, because there are often numerous hot gas-deficient 
mass substructures nearby to confuse matters and make de- 
termining which one the X-ray substructure originated from 
impossible (here we also have the time sequence to help us 
with this). 

Scenario 3 is sensitive to the choice of X-ray catalogue 
so, treating this type of detection as unwanted, we can con- 
clude that the Icr catalogue suffers from more false detec- 
tions (by our definition), which goes part way to explain its 
lower overall matching success. Scenario 2 can also be sensi- 
tive to the catalogue choice, as if the displacement between 
substructures is small, then the increase in area of a la X- 
ray detection can be enough to meet the overlap criterion in 
cases where it wasn't met for the 3cr X-ray detection. For this 
reason we continue to show the main results for matching 
to both the la and 3(t X-ray catalogues, although it should 
be noted this effect doesn't have a big impact on the match- 
ing success. Furthermore, in Section [71 where we simplify 
the discussion by showing results for one only X-ray cata- 
logue, it is the la that is chosen (despite its more numerous 
spurious detections) since it provides the most conservative 
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Figure 13. Matching success per area bin of 3cr 2D mass sub- 
structure catalogue to Icr (solid) and 3(t (dashed) 2D X-ray sub- 
structure catalogues versus area, for substructures above the pu- 
rity threshold (Agub = 10~^/i~'^Mpc^). Bins are equally spaced in 
log(ylgu]-|). Top: Fractional matching success when starting with 
the 2D total mass substructures and seeking a 2D X-ray counter- 
part for each. Bottom: Fractional matching success when starting 
with the 2D X-ray surface brightness substructures and seeking 
a 2D mass counterpart for each. 



estimate of the number of 2D mass substructures for which 
an X-ray counterpart cannot be found. 



6.2 Matching 2D X-ray substructures to 2D mass 
substructures with a 3D counterpart 

Subhalo mass is expected to be a crucial factor when look- 
ing for mismatches between DM and hot gas, as gas strip- 
ping procedures hav e more effect on low-mass subhaloes 
dTormen et al.l 120041 '). Here, the 3D subhalo data becomes 
an invaluable tool, not only because it effectively calibrates 
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our 2D total mass substructure catalogues, but because it al- 
lows us to probe the effect of subhalo mass on the success of 
the matching to the 2D X-ray substructures. We repeat the 
matching procedure outlined above, but this time only use 
2D mass substructures which have successfully been associ- 
ated with a 3D subhalo in Section[5l It should be noted that, 
in this section, 2D mass substructure catalogue now refers to 
the calibrated version of the original catalogue, containing 
only those 2D mass substructures which have a 3D subhalo 
counterpart. In the case of 2D mass substructures that have 
been associated with more than one subhalo (see Section [5]), 
Msub refers to the combined DM mass of these subhaloes. 

Table [2] shows the overall statistics for matching both 
the full and pure (i.e. ^sub.TM > 10""^ /i~^Mpc^) catalogues 
of 2D mass substructures with a 3D subhalo counterpart 
to the 2D X-ray substructure catalogues. It is worth high- 
lighting here that the results in Table [T] and [5] are almost 
identical for the pure 2D mass substructure catalogue be- 
cause, by definition, the vast majority of 2D mass substruc- 
tures in the pure sample have a 3D counterpart. Note that 
the ratios of X-ray substructures to total mass substruc- 
tures in the first row (full 2D mass catalogue) are larger 
than the equivalent ratios in Table [T] which does not in- 
clude any reference to the 3D subhalo data. This difference 
can be directly attributed the fact that a 3D subhalo coun- 
terpart could not be identified for around 10 per cent of the 
original 2D mass substructures (in the 3ct catalogue) and 

so A'^TM,Tabi<[2]/^TM,Tabi<[T] — O'^' whereas A^sb remains the 
same. 

There is a slight improvement in the fraction of 2D 
mass substructures matched to X-ray after calibrating 
the mass substructures against the 3D subhalo data (i.e. 

/matchcd,Tabi<l2] > /matchcd,Tabi<[lj) ' Primarily becausc this pro- 
cess will have removed any spurious detections from our 
2D mass catalogues. These are unlikely to be matched to 
an X-ray substructure, simply because they are due to dis- 
creteness noise in the total mass map or an artefact of the 
unsharp- masking procedure and, as such, we would not ex- 
pect these to be correlated with features in the X-ray surface 
brightness map. The removal of these unmatched mass sub- 
structures results in a boost to the overall matching success 
and slightly reduces the fraction of 2D mass substructures 
matched more than once (/muitimatchcd) to 3cr X-ray sub- 
structures. Despite this effect, however, the overall match- 
ing success still remains surprisingly low, with a maximum 
value of 45 per cent, achieved when matching to the Icr X- 
ray catalogue. 

This overall statistic is dominated by substructures with 
low associated 3D subhalo masses as these are far more nu- 
merous. From Fig. |3]we can estimate that there are approxi- 
mately 3 times more subhaloes with lO^'^ < Msuh / h~^}AQ < 
10^^ than subhaloes with Msub/Zi'^M© > 10^^ (for sub- 
haloes with their most bound particle within rsoo only). 
From the middle panel of Fig. [J] it is apparent that 85 — 95 
per cent of 3D subhaloes in this mass range (and with most 
bound particle within rsoo) have no hot gas. With these two 
results in mind, it is not surprising that the total percent- 
age of 2D mass substructures with 3D subhalo counterparts 
which also have X-ray counterparts is biased so low. This ef- 
fect can be further demonstrated by considering the overall 
matching success to X-ray for the pure 2D mass substruc- 
ture catalogue (second row. Table [Sjl . This is significantly 



Table 2. General results of matching between the 2D total 
mass substructures (which have been successfully matched to a 
3D subhalo) and 2D X-ray surface brightness substructure cat- 
alogues. Top rows: matching between all substructures in the to- 
tal mass catalogue, to all those in the X-ray catalogue. Bottom 
rows: matching between a pure (Agub xM > 10~^/i~^Mpc^) to- 
tal mass catalogue and all substructures in the X-ray catalogues. 
Columns: Multiple of cr used in indicated catalogue (X), Ratio of 
number of substructures ( ff^^ ), fraction of total substructures in 
the catalogue indicated that are successfully matched (/matched): 
fraction of matched substructures in the catalogue indicated that 
are matched more than once (/muitimatchcd)- 



^TM ^SB iV^TM fmatchcd fmultimatchcd 

Full 2D total mass substructure catalogue 

3 1 1.03 0.45 0.06 

3 3 0.62 0.34 0.07 

Pure 2D total mass substructure catalogue 

3 1 2.84 0.77 0.10 

3 3 1.72 0.67 0.10 



higher, with a maximum value of 77 per cent (again for 
the la X-ray catalogue). Here, on removing substructures 
with y4sub,TM < 10~^ /i~^Mpc^ to achieve a pure sample 
we have, by virtue of the Mgub — A^ub correlation (Equa- 
tion [12]), removed substructures with an associated value of 
A/sub < 10"/i-'Mq. 

It is clear from the overall matching statistics that the 
matching success depends heavily on the 3D subhalo mass, 
so we now examine this dependency in more detail. Fig. 1141 
shows the success in matching 2D total mass substructures 
to 2D X-ray substructures as a function of the associated 
subhalo mass. It should be remembered that Msub is the 
subhalo DM mass, not total mass, and therefore is inde- 
pendent of the amount of gas removal a subhalo may have 
undergone, other than the secondary effect of the remaining 
DM being more prone to tidal stripping. We have opted to 
include only those mass substructures matched to subhaloes 
above Msub = 10^ ^"'^M© (the completeness limit). How- 
ever, we note that the alternative choice of a sample with 
^sub,TM > 10~^ /i~^Mpc^ (the purity limit) made little dif- 
ference to this figure for the mass range shown. 



For the la X-ray catalogue, the matching success per 
mass bin rises gradually with DM subhalo mass: it is > 95 
per cent for cluster cores (Msub ^ lO^'*ft~^M0), ~ 95 per 
cent for groups (Msub ~ lQ^^h~^M.Q) and ~ 65 per cent 
for galaxies (Msub ~ lO^^/i~^M0). A similar trend is seen 
for the 3(7 catalogue, except that the success within a given 
mass bin is around 10-15 per cent lower. 

It is expected, based on the trend for /gas decreasing 
with subhalo mass already demonstrated, there would be 
a cut-off mass below which the fraction of substructures 
with an X-ray component would fall off sharply from 100 
per cent (similar to that seen in Fig |S] when matching 2D 
substructures to 3D subhaloes). This feature is present (at 
Msub ~ 10^^/i~^Mq), however the fall off is much more grad- 
ual and surprisingly the success rate is slightly below lOO per 
cent even above this value. The position of the cut-off mass 
and the rate of fall-off thereafter, is independent of X-ray 
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Figure 14. Fractional matching success per mass bin of 3(t 2D 
mass substructure catalogue to la (solid) and Scr (dashed) 2D 
X-ray substructure catalogues versus subhalo mass. Data are for 
2D mass substructures for which a 3D subhalo counterpart was 
found. Bins are equally space in log(Mi3ui,), the subhalo DM mass. 



catalogue choice and there is a deviation in the predicted 
success of only ~ 10 per cent. 

In order to link these results to the composition of the 
underlying 3D subhaloes, we examine the hot (T > lO^K) 
gas fractions, /gas, of the subhaloes that are matched to 
the 2D mass substructures. In cases where more than one 
subhalo is associated with the same 2D mass substructure, 
we calculate /gas for the most massive subhalo (but note that 
the exclusion of 2D substructures matched to more than one 
subhalo from the following analysis makes little difference to 
the results). Fig. llSl shows the average hot gas fraction, /gas, 
per mass bin for different samples of 2D mass substructures 
(binning is identical to that in Fig. I14|l . Interestingly, the 
average /gas for all 2D mass substructures (dot-dash) shows 
the same trend with subhalo mass as the matching success, 
suggesting that the two are closely linked, as it would have 
been reasonable to assume. 

The solid and dashed lines show the average /gas 
value for 2D mass substructures which are matched to a 
substructure in the Icr or 3cr X-ray catalogue, respectively. 
These values are slightly higher than those for all the 
substructures (with and without X-ray counterparts) 
suggesting that only low mass 3D subhaloes with a higher 
than average hot gas fraction will be successfully detected 
in 2D and matched in the X-ray maps. We also examined 
the minimum /gas of a total mass substructure with a 2D 
X-ray counterpart and the maximum /gas of one without a 
2D X-ray counterpart, per mass bin, for both X-ray cata- 
logues. These quantities effectively give the hot gas fraction 
thresholds that define the X-ray substructure catalogues. 
The minimum /gas was catalogue independent, suggesting 
the lower detection limit is dominated by another factor, 
however, the maximum /gas was found to be significantly 
higher for the 3a" catalogue at low masses (compared to 
the la catalogue and the average value), occasionally by 
a factor of around 2. Since the average /gas of a mass 
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Figure 15. Average hot gas fraction per mass bin of most massive 
3D subhalo matched to a 2D mass substructure in the 2D cata- 
logue (dot-dash) . Gas fractions are in units of Ob /Om • Other lines 
show the same quantity, but only for those 2D mass substructures 
that are also successfully matched to a 2D X-ray substructure. Bin 
values and other line-styles are as for Fig. 1141 



substructure successfully matched to an X-ray substructure 
is catalogue independent, yet the maximum /gas can be 
much higher for the 3ct catalogue, this suggests that that 
the difference in numbers of substructures contained in each 
catalogue is primarily due to hot gas distribution rather 
than mass. The objects picked up in the Icr catalogue but 
not in the 3a catalogue have /gas values higher than the 
average detected, i.e. if /gas is the controlling factor, they 
should appear in both catalogues. However, if a subhalo had 
a significant hot gas fraction, but this had been displaced 
from its centre or the dense core of the hot gas had been 
disrupted and so the peak in X-ray emission was not so 
bright, this could explain the same substructure being de- 
tected in the Icr but not the 3cr catalogues (remember the a 
values refer to cuts in the residual X-ray surface brightness) . 



6.3 Summary 

We have attempted to match every 2D mass map substruc- 
ture to a 2D substructure in the corresponding X-ray map 
and have shown that there are numerous occasions when this 
is not possible, highlighting differences between substructure 
in the hot gas and DM components. The frequency of match- 
ing failures clearly increases with decreasing subhalo mass: 
a few per cent of cluster cores (Mgub > 10^'^ h~^MQ), ~ 5 
per cent of groups (Afsub ~ 10^^h~^MQ) and ~ 35 per cent 
of galaxies (A/gub ~ 10'^^/i~^Mq) do not have X-ray coun- 
terparts. Interestingly, we also find that around a half of the 
X-ray substructures detected don't have counterparts in the 
mass maps. As more joint weak lensing and X-ray studies 
are undertaken, we predict more 'dark haloes' will be found, 
with these discoveries not restricted to rare, merger events 
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involving high mass subhaloes but wiU occur frequently on 
the galaxy-mass scale. 



7 DISCUSSION 

The benefit of performing a cosmological simulation is that 
it will best mimic the complicated processes taking place 
during the formation of real galaxy clusters. However, by 
the same token, it can then be difficult to untangle the in- 
fluence of one parameter or physical process on the conclu- 
sions, from that of another. In this section, we investigate 
the effects of the main selection parameters (Section I7.ip 
and the main model parameters (Section 17. 3|l on Fig. 1141 In 
addition, we define several broad categories for the fate of 
a mass substructure's hot gas component, as viewed in 2D. 
By exploring a Case Study from each category (Section l7.2p . 
we attempt to illustrate how the overall picture of the corre- 
spondence between the total mass and X-ray maps, shown 
in Fig. 1141 is built up. In Section 17.41 we make a prelimi- 
nary assessment of the potential impact of analysing maps 
with noise and observationally achievable resolution on our 
results. 



7.1 Selection parameters 

Section [6] dealt with all the cluster maps as one data set, 
however, in reality there are selection parameters which 
come in to play when observing clusters. The two most sig- 
nificant are redshift and dynamical state, the effects of which 
we investigate below. 



7.1.1 Variations with redshift 

The role of redshift is generally important when examining 
any class of astrophysical object, as it cannot be assumed 
that a population will not evolve significantly. In the case 
of galaxy clusters, our understanding of this factor is 
particularly significant if they are to become robust probes 
of the cosmological parameters. In this section we divide 
our maps into two samples, ^ z ^ 0.2 and 0.5 ^ 2 < 1.0 
(chosen to contain an equal number of snapshots; 11), and 
test if there is any difference in the results. 

Fig. [TS] shows the fractional success of matching 2D 
mass substructures, that have been associated with a 3D 
subhalo, to 2D X-ray substructures versus subhalo mass, 
split into low and high redshift bins. There is a trend for 
higher matching success between substructures in the total 
mass and X-ray surface brightness maps at higher redshifts 
for Maub < 10^^ft.~^MQ, most simply explained by the ar- 
gument that the subhaloes have had less time to suffer the 
effects of ram pressure stripping. 



7.1.2 Effects of dynamical state 

Since major mergers on the cluster mass scale are such en- 
ergetic events, it would be unsurprising if disturbed clusters 
exhibited more discrepancies in their hot gas and DM sub- 
structure. Indeed, some of the most extreme observational 
examples are found in highly disturbed c lusters, for exam- 
ple the bullet cluster (|Clowe et all \2004 ) and the ' cosmic 
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Figure 16. Matching success per mass bin of 3cr 2D mass sub- 
structure catalogue to la 2D X-ray substructure catalogues versus 
subhalo DM mass. Data divided by redshift: z 0.2 (dashed) 
and 0.5 ^ 2 < 1.0 (dot-dashed) Solid line shows combined data 
from Fig. 1141 



train- wreck' in AbeU 520 (|Mahdavi et al.ll2007l ). There is 
also much debate about the significance of the effect that 
merger activity has on bulk properties of galaxy clusters 
potentially making it an important selection effect. Simula- 
tions of isolated cluster mergers have suggested that massive, 
correlated luminosity and temperature boosts are associated 
with major mergers and that these could cause the masses 
of high-redshift clusters to be overe stimated from both the 
M - Tx and the M - Lx relation (jRicker fc SarazinllioOll : 
iRandall et"ai1l2002l '). It has also been suggested that, if this 
effect is real, such systems would stand out from scaling 
relation s (lO'Hara et al."2006'). but in recent high resolution 
studies l|Poole et al.li200 6. 2 0071 . 1200^ and cosmological sim- 
ulations where multiple mergers occur, no such simple cor- 
relation between scatter in the Lx — Tx relation and visible 
evid ence of ongoing ma j or merger activ ity has been found 
(e.g. iRowlev et ai]|2004l : iKav et al.ll2007l ). 

In order to divide our images into just two subsets (ma- 
jor merger or not), it is necessary to have an additional 
technique to calibrate the centroid shift variance (described 
in Section 12. 3p and determine which value of this statis- 
tic marks the threshold between these two states. Since we 
are interested in separating out the most extreme merging 
events, as this should make any trend stand out, we choose 
a value of the centroid shift variance which singles out the 
highest peaks in this quantity, which we determine to be 
0.1. The sample is then split into two - those snapshots with 
values above the threshold and those with values below and 
this division is confirmed by examination of X-ray surface 
brightness contour maps (see Fig. [T]). 

Fig. [17] shows the matching success for the disturbed 
(dashed) and relaxed (dot-dash) samples, versus subhalo 
mass. It is clearly evident that all the matching failures 
above Msub — 3 x 10^'' lie in the disturbed sample, and visual 
inspection of the maps confirms they are all undergoing sig- 



20 Leila C. Powell et al. 



" 0.8 - 



^ 0.6 -_. 



D 

E 

o 
c 
o 

tj 
o 



0.4 



0.2 



0.0 



2.0 2.5 3.0 3.5 



4.0 



4.5 

)] 



5.0 



Figure 17. Matching success per mass 6m of 3cr 2D mass sub- 
structure catalogue to la 2D X-ray substructure catalogues ver- 
sus subhalo DM mass. Data split into two subsets according to 
dynamical state (see text for details): relaxed (dot-dashed) and 
disturbed (dashed). Solid line shows combined data from Fig. 1141 



nificant mergers or collisions. The trend is actually reversed 
for low masses; substructures here have a lower probability 
of having an X-ray counterpart in relaxed clusters. This can 
be explained in the same way as the dependency on red- 
shift in this mass range. Substructures in relaxed clusters 
have been there longer since, by definition, the last merger 
event was some time ago and therefore have been subject 
to stripping processes for longer. For comparison with our 
two redshift samples, we compute the mean redshift of the 
snapshots in our disturbed and relaxed samples, which we 
find to be 0.55 and 0.37 respectively. This highlights that 
there is some degeneracy between the effects of redshift and 
dynamical state, although these values are much closer than 
the mean redshifts of our redshift samples (0.72, 0.10; snap- 
shots are equally spaced in time not redshift) and so the 
fact we get such a dramatic difference suggests looking at 
redshift alone is not sufficient. 



7.2 Detachment of hot gas from dark matter 
subhaloes 

In the previous section we have investigated the overall 
probability of finding a 2D X-ray counterpart for the 2D 
total mass substructures we detected. In order to fully probe 
all the factors which result in a matching failure in our 
analysis, a more detailed treatment (e.g. careful tracking of 
subhaloes between snapshots with higher time resolution 
or idealised simulations of individual mergers) would be re- 
quired. Nevertheless, it is very informative to examine some 
of these matching failures in more detail to gain insight 
into the variety of scenarios that occur. It is reasonable to 
assume that ram pressure stripping is the main culprit in 
the removal of hot gas, however it is interesting to note that 
there are several distinct realisations of the outcome of this 
process in the maps. To illustrate these, we choose a small 



subset of substructures by visual inspection and follow 
these up in 2D and 3D. We present these case studies below. 



7.2.1 Evidence of partial ram pressure stripping. 

Although we are primarily concerned with cases in which 
matching between the 2D total mass and 2D X-ray surface 
brightness substructures fails, it is interesting to note that 
we observe the signatures of ram pressure stripping in 
objects for which the match is still achieved. 

Case Study 1. Fig. [18] shows the development of a 
large tail of hot gas which streams behind the main sub- 
structure as it passes near to the cluster centre (its 3D 
physical displacement in the middle panel is approximately 
0.5r5oo). The substructure is indicated with an arrow and 
has Msub ~ 1 - 3 X lO^^/i"^M0. Note the temporary de- 
crease in mass (about a factor of 2) in the middle snapshot. 
This is a well-known issue with SUBFIND, where subhaloes 
become more difficult to distinguish when in close proxim- 
ity to the cluster's c entre, due to their low density contrast 
(jLudlow et al.ll2009D . In this case, around half of the sub- 
halo's DM particles in the left-hand image are deemed to 
belong to the main cluster in the middle image, decreas- 
ing to around 20 per cent in the right-hand image, as the 
subhalo moves away from the core region. 

The substructure is clearly visible in both the X-ray 
image and the mass contours and is detected and matched 
in 2D at each time shown, regardless of choice of the la 
or 3(T X-ray catalogue. The SUBFIND data reveals that the 
mass of gas associated with the corresponding subhalo is 
reduced to 60 per cent of its initial value over the course of 
this time sequence, however. The subhalo may well be on 
a highly elongated orbit and may eventually be depleted of 
enough of its hot gas such that it is no longer detected in 
the X-ray surface brightness map. As such, it is likely that 
this represents what is, for many substructures, the first 
stage in a time sequence that eventually leads to matching 
failure. However, this stage may not be visible for many 
substructures, depending on the inclination of the gas tail 
with respect to the line of sight and the amount of X-ray 
emission from the stripped gas. 

This case study also demonstrates that our detection 
technique is not too sensitive to the stripping of the outer 
regions of a substructure's hot gas and that matching fail- 
ures therefore represent an extreme depletion or complete 
removal (or displacement) of the X-ray emitting component. 



7.2.2 Matching failure due to complete ram pressure 
stripping. 

We now examine the scenario whereby stripping of a 
substructure's hot gas results in the DM component turning 
up as an unmatched total mass substructure. There are two 
distinct categories here: matching failure due to the spatial 
displacement of the hot gas component and matching failure 
due to the hot gas component being erased completely from 
the X-ray surface brightness map. 
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Figure 18. Sequence showing strong evidence of stripping in cluster B for 0.16 $5 ^ ^ 0.21, but neither detection nor matching procedures 
fail. Image is composed of the logarithmic X-ray surface brightness map with X-ray surface brightness (white) and surface mass density 
(black) contours, equally spaced in log, overlaid. Inner contours have half the spacing of outer ones in order to highlight structure in the 
core region. M^ub ^ 2.5 X 1.3 X IG^/i-IMq and 2.5 X lO^^h'^MQ, respectively 




Figure 19. Sequence showing progressive separation of 2D total mass substructure and its 2D X-ray counterpart for cluster B at 
0.05 < z < 0.01, which results in a matching failure in the final map of the sequence. Image is composed of the logarithmic X-ray 
surface brightness map with X-ray surface brightness (white) and surface mass density (black) contours, equally spaced in log, overlaid. 
A/sub 2.6 X lO^2/i~^M0,2.4 x lO^^h-^M© and 3.0 X lO^^/i^^M©, respectively 



Case study 2. |ProkIiorov fc Purred (|2007l ) produce an 
analytical model describing the increasing separation be- 
tween the DM and hot gas components of a substructure 
moving through the ICM. Indeed this type of displacement, 
where the hot gas substructure remains intact yet is clearly 
displaced from the DM, is another reason for matching fail- 
ures in our analysis. Fig. [19] shows a sequence of maps in 
which an X-ray counterpart is initially found for a total mass 
substructure, but then clearly becomes spatially separated 
to the point where matching fails. It is clear that should 
this type of displacement occur along the line of sight, we 
would not be aware of it, however this isn't a shortfall of our 
method; we want to be subject to the same restrictions as 
observers. In the left panel the substructure is well matched 
to its X-ray counterpart (indicated by one arrow only as 
the contours and surface brightness peak are coincident). In 
the middle panel, due to the effects of ram pressure, the 
X-ray component is slightly offset, yet matching is still suc- 
cessful for both the 3a and la X-ray catalogues. Finally, in 
the right panel, the progression of the X-ray component has 
been slowed so much that it is significantly displaced from 



the total mass substructure and is, therefore, not matched 
regardless of the choice of X-ray catalogue. 

A 3D analysis of this time sequence confirms this pic- 
ture. The bound gas mass of the corresponding subhalo de- 
creases to almost zero between the left-hand and middle 
panels of Fig. [19] and the subhalo is completely gas-free by 
the right-hand panel. We identify, in the right-hand panel, 
the location of the gas particles which originally belonged 
to the DM subhalo and confirm that these exist as a clump 
which is coincident with the substructure in the X-ray sur- 
face brightness image. 

This separate hot gas component will eventually be 
completely disrupted, leaving the mass substructure (which 
will remain intact longer) with no trace of an X-ray counter- 
part. This situation is seen in the maps fairly frequently and 
presumably the separation procedure described above is the 
precursor to this. However, we also encounter an example 
of more immediate erasing of substructure from the X-ray 
surface brightness map, which we describe in the next case 
study. 

It should be noted that the survival time of the 
gaseous component of subhaloes has been shown to be 
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dependent on the numerical techniques employed and the 
resulting success with which hydrodynamica l instabilities 
that expedite gas stripping are captured. lAgertz et alJ 
|20o3) show that SPH (with standard artificial viscosity) 
cannot capture Kelvin-Helmholtz Instabilities ( KHI) as well 
as Ad aptive Mesh Refinement codes. Indeed, iDolag et al.l 
l|2008l ) demonstrate that using a low-viscosity scheme 
(less damping of the KHI) in GADGET-2 results in smaller 
gas fractions for subhaloes inside a cluster's virial radius, 
suggesting such issues will impact on studies such as this 
one. More recently, improvements to the SPH methodology 
have been suggested, which increase its ability to cap ture 
KHI (|Pricdl2008l : iKawata et ai]|2009l : [Read et aLllioog ). It 
is therefore a matter for further investigation, whether the 
concentrations of stripped subhalo gas which give rise to 
X-ray surface brightness substructures, as illustrated in this 
case study, would still occur when the growth of KHI are 
simulated reliably. 

Case study 3. Fig. 1201 shows the encounter of a sub- 
structure (Afsub — 10^'^/i~^Mq) with the inner regions of 
the main cluster at z « 0. In the left panel the substruc- 
ture is detected (and clearly visible by eye, indicated with 
an arrow) in both the total mass and X-ray surface bright- 
ness maps and the matching procedure is successful. In the 
middle panel there is a clear double-peaked structure in the 
centre of the X-ray map, yet this is absent in the total mass 
contours and the substructure is not detected as a separate 
object from the cluster core in 2D. SUBFIND individually re- 
solves the subhalo in 3D, albeit with a significantly reduced 
mass due to the effect discussed in Case Study 1. Therefore, 
although we have two 3D haloes (the main cluster and the 
subhalo), their proximity means both are attributed to the 
same 2D mass substructure; essentially the resolution limit 
of our detection procedure has been exceeded here. 

In the right-hand panel the substructure is again de- 
tected in the total mass map (its position in the contours 
is indicated with an arrow) but there is no corresponding 
detection in this region of the X-ray map. The DM sub- 
structure has been completely stripped of its hot gas and is 
also stripped (tidally) itself; according to SUBFIND the cor- 
responding subhalo has no bound gas particles in the final 
panel and the DM mass has been reduced to roughly 40 per 
cent of its value in the first panel. Indeed the rest of this 
subhalo's DM and gas particles are found to belong to the 
main cluster at the end of this time sequence. 

Note, however, the edge-like feature present instead 
which appears to lag behind the 2D mass substructure (de- 
tections of this are made in both X-ray catalogues but are 
not matched to the mass substructure as the positions are 
not coincident). A tail of stripped hot gas is also visible on 
the right-hand side of this map which resulted from the sub- 
structure's approach and is also apparent in the left and mid- 
dle panels. A detailed discussion of the gaseous f eatures that 
can a rise during such interactions appears in IPoole et al.l 
(|2006l ). 

We also note that there are similar scenarios whereby 
a collision results in the disruption of the X-ray emitting 
gas, rather than its complete removal as seen here. In these 
cases, while there is still evidence of gas in the vicinity of the 
mass map substructure, a defined peak is no longer visible 
and the detection of an X-ray substructure can then be cat- 



alogue dependent (difi'erent catalogues impose different cuts 
on the residual surface brightness). The issue of catalogue 
dependence is discussed in the following case study. 

If we examine the maps immediately preceding this 
sequence we observe the same substructure undergoing 
stripping of its hot gas on an earlier passage through the 
cluster's inner regions, adding weight to our supposition 
that Case Study 1 may be the precursor to the hot gas 
being removed completely. In fact, it is plausible that all 
the case studies may simply be different moments in a 
sequence which substructures that continue to orbit within 
the cluster long enough are subject to: stripping of outer 
regions of hot gas, displacement of remaining hot gas then 
complete disruption of the hot gas substructure. 

Case study 4. We now examine a scenario in which a 
close encounter between the main cluster and a subhalo has 
resulted in the hot gas from one object being removed and 
assimilated into the main cluster's ICM. However, in this 
case the conglomerate hot gas does not correspond spatially 
with a mass map substructure. This case study highlights 
how challenging it can be to correctly determine which X- 
ray substructures are associated with which mass substruc- 
tures when working in 2D and the limitations of our current 
detection and matching schemes. 

In Fig. [21] we present an image of cluster A at 2: = 1 
(middle panel) which bears similarity, in terms of i ts con - 
figuration, to the recent observation in lBradaF et all l|2008h . 
This observation is on a larger scale, however, depicting the 
merger of two clusters, with nearly equal masses. The left 
and right-hand panels show the maps for the snapshot di- 
rectly before and after, respectively, to provide some insight 
however, due to the complexity of the merger, a detailed 
subhalo merger tree would be required to unravel the full 
series of events. Here we focus on the middle panel. 

It is clear that the inner region consists of three DM 
haloes (all detected in 2D by our algorithm) and only two X- 
ray peaks (also detected). The western X-ray peak is coinci- 
dent with (and successfully matched to) one of the DM sub- 
structures, however, the eastern X-ray peak lies in-between 
the other two substructures. The smaller of these DM sub- 
structures (Maub — 1.6x lO^'^/i~^M0) is matched to the east- 
ern X-ray substructure when using either the la or 3a" X-ray 
catalogue. We can see, however, that this substructure is still 
offset from the bulk of the hot gas in this region and instead 
just overlaps with a tail of material extending outwards to- 
wards it. Indeed, in 3D this subhalo has a very low gas frac- 
tion, suggesting its match with the X-ray substructure is pri- 
marily a projection effect. The larger of these DM substruc- 
tures appears completely devoid of hot gas and is the main 
cluster as defined by SUBFIND [Ms^h — 1-9 x lO^*/i~^M0). 
The displacement of this surface mass density peak from 
the eastern X-ray surface brightness peak is 90/i~^ kpc. It is 
interesting that while the projected mass peak of the main 
cluster is offset from the bulk of the X-ray emission, its most 
bound particle actually coincides with the maximum X-ray 
surface brightness (the white spot in the image). When we 
use the la X-ray catalogue, the main cluster core is actually 
matched to a very small, separate X-ray substructure and 
so does not show up as an unmatched object in this data 
set. This pair would add to the scatter in the area-area cor- 
relation plot (Fig. I12p suggesting that, in future work, poor 
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Figure 20. Sequence stripping of X-ray gas from a substructure in for cluster A for 0.07 ^ ^ ^ 0., matching fails in the last image 
of the sequence. Image and contours as in Fig. [T9l A4ub - 1-0 X lO^^/i-^M©, undefined (see text for details) and 4.4 X IO^^/j-IMq, 
respectively. 




Figur e 21. Image of cluster A undergoing a merger at z ~ 1, with the middle panel bearing remarkable similarity to that in lBradac et al.l 
1120081 ). Left-hand panel: Masses of three main inner substructures are 2.8 X 10^^/i~^Mq (smaller peak), 1.8 X IO^^/x^^Mq (large eastern 
peak) and 3.7 X IO^^/i^^Mq (large western peak). Middle panel: The western X-ray peak is coincident with one of the DM substructures 
yet the eastern X-ray peak lies in-between the other two substructures. The larger of these appears completely devoid of hot gas and is 
the main cluster as defined by SUBFIND {M^ub — 1-9 x lO^'*/i~'^M0). The displacement of this surface mass density peak from the eastern 
X-ray surface brightness peak is 90h~^ kpc. The smaller DM substructure has M^ub — 1-6 X lO^'^/i^^M©. Right-hand panel: Masses of 
two main inner substructures are 2.0 x l0^^h~^yiQ (large eastern peak) and 2.1 x 10^*/i~^Mq (large western peak). Image and contours 
as in Fig. 1181 Note that north is up and east is left. See text for further details. 



agreement in area of matched pairs could be used to remove 
dubious matches. 

This case study illustrates that the use of the 3a X- 
ray catalogue is the most likely to allow retrieval of all high 
mass substructures with no significant hot gas component 
as these will fail to be matched. The la X-ray catalogue is 
less stringent and will detect smaller amounts of gas and will 
also result in the matching of substructures that are slightly 
offset, since the X-ray substructures in this catalogue are 
more extended. However, by the same token, it is also more 
prone to false matches than the 3a" catalogue. 



7.3 Gas physics model parameters 

Our main results focussed on a set of non-radiative clus- 
ters, the simplest model for the ICM within a cosmological 
context. It is well known, however, that additional physical 
processes must operate in clusters; scaling relations such as 
the X-ray luminosity-temperature relation are different to 



what is exp ected from the so-called self-similar model (e.g. 
lKaiseij[l99ll ) . The most favoured explanation for the altered 
similarity of clusters is that the ICM has undergone an in- 
tense, and perhaps extended, period of heating due to galac- 
tic outflows (from stars and active galactic nuclei) . Radiative 
cooling also plays a role, selectively removing the low en- 
tropy gas, although is completely reliant o n subsequent heat- 
ing to avoid a cooling cata strophe (e.g. iBabul et al. I I2OO2I : 
iMcCarthv et al.l |2004 120081 ) . 

Investigating the full effects of cooling and heating is 
beyond the scope of this paper, but we have performed a 
preliminary investigation on the effects of cooling on our re- 
sults. To avoid over-cooling the gas and motivated by the 
observations that stellar populations in clusters are old (e.g. 
Thomas ct al. 2005), we only allow the gas to cool radia- 
tively at early times, until a reasonable fraction of gas has 
cooled and fo rmed stars . We a dopted the same procedure 
as outlined in iKav et all l|2004 ). assuming a metal-free gas 
(Z — 0). The simulation of cluster A was repeated and 
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Table 3. Properties of clusters A (NR and CL denote the non- 
radiative and cooling runs respectively), B and C at 2; = 0. M500, 
Tgi and Lx.sim a^re calculated from the simulation data directly. 
Lx,obs a^rs the observed X-ray luminosities for clusters with the 
same value of M500 or T^i, calcu lated using the L x — A^soo or 
the Lx — Tx scaling relation from lPratt et al" ] l|2009h . denoted by 
L-M and L-T, respectively. 



Cluster 
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/^X,sim 
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22.3 
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cooling was switched on until « ~ 5, when around 10 per 
cent of the gas in the high-resolution region had formed 
stars, in agreement with near-infrared observations, which 
suggest that the stellar mass, as a percentage of the ICM 
gas mass, is around 10 per cent on avera ge (e.g. iLin et all 
I2OO3I : iBalogh et al.ll200l[|Cole et al.ll200lh . The cooled frac- 
tion in the cluster, within rsoo, at 2 = is about 2 per cent, 
slightly higher than observed bv lLin et al.l l|2003l ). 




log,o[Msub(h"' 10^° Me)] 

Figure 22. Cumulative subhalo DM mass functions for cluster A 
cooling run (solid line), and cluster A non-radiative run (dashed 
line), at 2 = 0. Data are for subhaloes with their most bound 
particle within rsoo- 



We use obse rved Lx — M500 and Lx — Lx relations 
l|Pratt et al.ll2009l ') in order to compare the simulated X-ray 
luminosity (Lx,sim) of cluster A at 2 = with that expected 
from observations (Lx.obs) based on both its M500 value 
and i ts spectroscopic-like temperature, Tsi (jMazzotta et al.l 
I2OO4I ). Table [3] summarises these properties, including the 
values for the other clusters as a point of comparison. The 
non-radiative version of cluster A has a similar X-ray lu- 
minosity to that observed based on its mass, but a much 
higher luminosity than observed (Lx,sim/Lx,obs ~ 5) based 
on Tsi. Similarly, clusters B and C also exhibit X-ray lumi- 
nosities close to those observed for their masses at z = 0, 
yet are very over-luminous for their temperatures. Assum- 
ing the observed mass determinations are accurate, then the 
main difference is that the simulated Tgi is too low in the 
non-radi ative model due to the presence of too much cold 
gas (see lKav et"ar '2008Y Turning on high-redshift cooling 
in cluster A preferentially removes this cool gas, bringing 
down the luminosity (Lx.sim/Lx,obs ~ 0.6, based on mass) 
but increasing Tsi such that Lx,sim/Lx.obs ~ 1-6, based on 
temperature. So, overall, the cooling model is closer to the 
observed (best-fit) Lx — Tx — M500 plane. In future, we will 
also consider the additional effects of heating from super- 
novae and active galactic nuclei. It will be interesting to see 
how the competing effects of cooling and heating affect the 
structure of the subhaloes. 

The cluster identification procedure is the same as de- 
scribed in Section [2j except that in order to ensure we follow 
the same object in all resimulations of the same cluster, the 
list of cluster candidates in the cooling run is searched for 
the best match to the selected object in the non-radiative 
run. In Fig.[2]the surface mass density (left) and X-ray sur- 
face brightness (right) maps for the cooling run (first row) 
and non-radiative run (second row) of cluster A at 2; = 
can be compared. Qualitatively, the two sets of maps ap- 
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Figure 23. Average hot gas fraction per mass bin (in units of 
f^b/f^m = 0.15) for cooling run (solid) and non-radiative run 
(dashed) of cluster A. 



pear similar, suggesting that cooling at high redshift does 
not strongly affect our main results. 

We can once again examine the properties of our 3D 
subhalo sample, to see the role cooling has played. Fig. 1221 
shows the 2 = subhalo DM mass function for the cooling 
(solid line) and non-radiative (dashed line) runs of cluster A. 
The DM mass functions agree well, although it is apparent 
there are more low DM mass subhaloes (a similar effect is 
also seen at 2 = 0.5 and 1), suggesting that the central 
condensation of the baryons deepen the potential wells and 
reduce the amount of disruption of the DM. 
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We examine the effect of cooling on the hot gas within 
the 3D subhaloes, by computing the average hot (T> 10*' K) 
gas fractions within each mass bin. The results are shown 
in Fig 1231 for the non-radiative (dashed) and cooling (solid) 
versions of cluster A (expressed in units of the global value, 
f2b/fim = 0.15). In the non-radiative cluster, the hot gas 
fraction increases with subhalo mass, reflecting the increas- 
ing ability of the subhaloes to retain their hot gas as their 
potential wells deepen (note the main haloes are also shown). 
In the cooling run, the same trend is seen, but the gas frac- 
tion is lower at all masses (compared to the non-radiative 
run). This is due to cooling causing the additional depletion 
of the hot gas reservoir by transforming it into cold gas (and 
eventually, stars). 

As a result of their shallower potential, ram pressure 
stripping of hot gas is most effective in the low mass (~ 
10^^ /i~^Mq) subhaloes, as indicated by their very low av- 
erage gas fractions (~ 0.2) in both runs. It is at this mass 
range, however, that cooling is also most effective because 
these objects form at high redshift and have short cool- 
ing times. We see the result of this effect when comparing 
the average subhalo total baryon fractions in the two runs; 
the cooling cluster has around 40 per cent more baryons at 
Afsub ~ lO^^/i~^M0, yet the total baryon fractions agree 
well between the runs at higher masses. 

The procedure to detect 2D substructures in maps of X- 
ray surface brightness and total mass density (described in 
Section[4| is applied, with identical parameters, to the cool- 
ing run of cluster A. As could be expected from the discus- 
sion of the 3D subhalo data for this run, more 2D total mass 
substructures are found; the 3a 2D mass catalogue contains 
611 substructures compared with 473 in the non-radiative 
run of the same cluster. The number of substructures in the 
la X-ray catalogue is very similar, with 395 compared to 
363 in the non-radiative cluster. 

The matching of 2D total mass substructures to the 3D 
subhalo data was undertaken for the cooling run as described 
in Section [51 and an assessment of completeness and purity 
of the 2D catalogue was repeated. The choice of the 3a 2D 
total mass catalogue was again deemed the most appropriate 
and the purity and completeness limits already established 
were found to be valid for application to the cooling run 
data. 

We now consider the likelihood of finding an X-ray 
counterpart for the substructures in the calibrated 2D total 
mass catalogue (i.e. the catalogue containing only those sub- 
structures that were successfully identified with a 3D sub- 
halo). Fig. 1241 shows the matching success to X-ray for total 
mass substructures. Surprisingly, the results for the cool- 
ing run (solid line) match those for the non-radiative run 
(dashed line) very closely. In particular, although Fig. 1231 
shows subhaloes of all masses are more depleted of hot 
gas in the cooling run, this doesn't seem to translate into 
a significant decrease in the likelihood of finding a 2D 
X-ray counterpart for 2D mass substructures, except for 
Msub ~ 1O^^^"^M0. OveraU, it seems that the introduc- 
tion of high-redshift cooling does not affect the main results 
and therefore that our non-radiative results are not too sen- 
sitive to the gas physics model employed (although further 
investigation into the effects of cooling plus feedback would 
be desirable). 
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Figure 24. Matching success per mass bin of 3cr 2D mass sub- 
structure catalogue to la 2D X-ray substructure catalogues for 
cooling run (solid) and non-radiative run (dashed) versus subhalo 
mass. 



7.4 Towards realistic observations 

The main results of this paper have focussed on 'perfect' ob- 
servations. Although a detailed analysis of all the potential 
observational and instrumental effects is beyond the scope 
of this paper, we now consider the impact on our results of 
using an observationally achievable map resolution and in- 
cluding basic noise in the maps. We defer a more detailed 
treatment of noise and instrumental effects to future work. 
This analysis is undertaken on the non-radiative simulations 
as this provides us with a larger sample of clusters and we 
have shown that the impact of high-redshift cooling is min- 
imal. 



7.4-1 Introducing noise and degrading the resolution. 

For the purposes of adding noise and adopting a realis- 
tic resolution, we opt to place the clusters at z = 0.2 as 
this is both the redshift at which our X-ray map resolu- 
tion can be achieved by XMM and is also a redshift rep- 
resentativ e of recent o bservations (the bullet cluster is at 
z » 0.3 JClowe et al.l [2004) and AbeU 520 is at z ^ 0.2 
iMahdavi et al.ll200'i1 )). 

As described in Section l4Tl maps of both types were first 
smoothed with a Gaussian kernel with FWHM — 15/i~^kpc, 
equal to the spatial resolution of the simulations. While this 
resolution is potentially achievable in X-ray observations, it 
is necessary for us to increase the smoothing slightly (to 
FWHM = 25/i~^kpc) in the presence of noise, but note that 
this resolution is still much higher than currently achiev- 
able with weak lensing analyses. In this case, the reso- 
lution that can be obtained is dependent on the number 
density of background galaxies: for ground-based data, this 
angular resolution is typically 1 arcminute, yet for space- 
base d data this can be im proved to around 45 arcseconds 
(see iHevmans et al.l |2008| . for example) . To investigate the 
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Figure 25. Images of cluster A at 2 0. First row: Pro- 
jected mass map witli resolution of 15h~^\qpc (left) and mass 
map with Gaussian noise and resolution degraded 100/i~^Kpc 
(right). Second row: X-ray surface brightness image with res- 
olution of 15/i~^kpc (left) and X-ray surface brightness image 
with 25/i~^kpc resolution and Poisson noise (right). Note that 
the noisy X-ray image would typically be more heavily smoothed 
for presentation purposes. See text for details of the noise models. 
This image is featured in Case Study 3, FiE l20l 



impact of this decreased resolution, we now adopt a pre- 
liminary Gaussian kernel with FWHM= 100/i~^kpc when 
analysing the projected mass maps (corresponding to ap- 
proximately 45" angular resolution at z = 0.2) and increase 
CT2 accordingly. 

We add Poisson noise to the X-ray maps by making 
the crude approximation that the photon number is pro- 
portional to the X-ray surface brightness. We find that 10^ 
photons (corresponding roughly to an exposure of 20 ks) al- 
lows us to recover the majority of substructures that were 
detected in the absence of noise. We add Gaussian noise to 
the mass maps with zer o mean and with a variance deter- 
mined bv Ivan Waerbekd l|2000l ). The latter is given by, 



47rCT^n£ 



2(T^ V 2 (T ^2a 



(15) 



for a pixel of size a in a weak lensing mass reconstruc- 
tion and is due to the intrinsic ellipticities (with rms, a^) 
of background galaxies with an average density of rig. Ecrit, 
the critical surface density for lensing to occur, is given by, 



47rGDoi,A. 



(16) 



where Doi is the angular diameter distance between the 
observer and the lens, Dos is that between the observer and 



the galaxies and Dis is that between the lens and the galax- 
ies. Note that we fix the redshift of the lens to be 2 = 0.2, as 
outlined above. We use typical values of rig = 100 galax- 
ies arcmin~^ (for space-based data) and o"e = 0.3 (e.g. 
IStarck et alll2006l ) and, we assume that the galajcy elliptic- 
ities were smoothed with a Gauss ian of standard deviation 
a = a prior to reconstruction, as in lPuchwein fc Bartelmamil 
( 20071 ). We note that the pre-smooth damps the noise a lit- 
tle, but our aim in this preliminary investigation into the 
impact of noise is simply to make an estimate of the noise 
level. 

Fig. illustrates how the original maps (left-hand col- 
umn) are afi'ected by the smoothing and addition of noise 
(right-hand column). While some small substructures are 
still visible in the X-ray map (albeit made less distinct by 
the noise), all but the largest substructure has been erased 
from the mass map. 

7.4-S The impact of noise and resolution. 

We now review our main findings in order to make a pre- 
liminary assessment of how they are affected by noise and 
degraded map resolution. 

First, we re-evaluate the relationship between the 2D 
mass map substructures and the underlying distribution of 
3D subhaloes. The impact of just degrading the map reso- 
lution is significant and reduces the number of subhaloes 
detected above 10^^/i~^Mq to around 60 per cent of its 
value in the original maps. When noise is also included, 
there are further detection failures, most frequently below 
5 X lO^^/i~^M0, which reduce the total number of detections 
above 10^^ h~^MQ by an additional 5 per cent. 

Fig. compares the completeness of the 3a 2D mass 
substructure catalogue obtained from the degraded resolu- 
tion, noisy maps (solid line) with that obtained from our 
high resolution, noise-free maps (dotted line, taken from 
Fig. IS}. We can see that the mass threshold for 90 per cent 
completeness (per mass bin) is now around an order of mag- 
nitude higher. For the 3cr catalogue, 90 per cent complete- 
ness is now achieved only above 10^^ h~^MQ and a factor 
~ 8 fewer substructures are detected in total (159 cf. 1233). 

Despite the impact of noise, a correlation between the 
area of the 2D substructure, A^uh, and the DM mass of its 
3D subhalo counterpart, Msub, is still evident in Fig. 1271 The 
reduction in the number of detections (immediately appar- 
ent in Fig. I27|l translates into higher la errors on the best- 
fitting line. The purity of the 3a substructure catalogue is 
very high (only 1 per cent are false detections), so we do 
not define a purity threshold here, but simply include all 
substructures in the fit. 

The normalisation and slope are now 4.7±0.2 and 0.9± 
0.1, respectively. The change in the former is most significant 
and can be attributed to the use of a larger kernel, resulting 
in a given 2D substructure having a larger area than before. 

We also re-examine the likelihood of finding an X-ray 
counterpart for all 2D mass map substructures in the 3a 
catalogue (the results for our fiducial data set are presented 
in Fig. I14|) . Above lO^^/i~^M0 (our completeness limit) we 
now find X-ray counterparts in the la catalogue for all of the 
2D mass map substructures. Using the 3a X-ray catalogue, 
however, we fail to find matches for a few high mass mass 
map substructures. The difference between the two X-ray 
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Figure 26. Fractional matcliing success per mass bin of 3D sub- 
lialoes to 2D substructures in the noisy, degraded resolution mass 
maps (solid) and the original mass maps (dotted) as a function of 
subhalo DM mass for 3cr 2D catalogue. Bins are equally spaced 
in log(A4ub)- 



matching failures typically occurred in complex merging 
cores which now, in some cases, cannot be individually re- 
solved. This can result in two merging cores being detected 
as one extended mass substructure, facilitating matching 
with an X-ray substructure. This is an issue that requires 
further investigation. The current criterion for a match is 
any degree of overlap between the mass and X-ray substruc- 
ture. Since substructures in the realistic mass map substruc- 
tures have a much greater spatial extent due to the lower 
resolution, they could be associated with an X-ray substruc- 
ture which is significantly offset from the mass peak. A more 
detailed follow-up of X-ray-mass matches in the context of 
more realistic maps would be an interesting extension to 
current work. 

While we have shown that only a few discrepancies be- 
tween substructure in the X-ray and the mass maps could 
be observed currently, our fiducial results show there is an 
abundance of these to be uncovered. A detailed substruc- 
ture comparison, such as the one undertaken here, will yield 
a wealth of interesting results when predicted improvements 
in lensing mass map resolution are achie ved. For example, 
a resolution of ~ 10/i^^kpc is forecast bv lCod (|200a ) based 
on a novel strong lensing analysis technique. 




-2.5 -2.0 -1.5 -1.0 
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Figure 27. Correlation between physical area, Agub, of 2D sub- 
structure in noisy, degraded resolution mass map (Sc catalogue) 
and DM mass of 3D subhalo, M^^^^, to which it is matched. Filled 
squares show SUBFIND background haloes (see text for details). See 
text for slope and normalisation of best-fitting line. 

catalogues here arises in situations where there is hot gas in 
the vicinity of the mass substructure, but it has been dis- 
rupted and so does not have a defined peak; the la catalogue 
detects this whereas the 3a does not. 

The other scenarios in which we no longer find matching 
failures (but did previously) is an effect of the reduced mass 
map resolution. At low subhalo masses we are now unable 
to resolve the 2D mass map substructure. At high masses. 



7.5 Summary of Discussion 

In this Section, we have discussed in detail the reasons for 
failing to find an X-ray counterpart for all of our 2D mass 
substructures. We have demonstrated two distinct scenar- 
ios that give rise to a 2D total mass substructure not be- 
ing matched to a 2D X-ray surface brightness substructure: 
spatial separation of the X-ray component and destruction 
(or disruption) of the X-ray component. We have also high- 
lighted the dependence of the matching procedure on choice 
of X ray catalogue. 

We have examined how several factors affect the likeli- 
hood of finding X-ray counterparts for substructures in the 
total mass maps. The inclusion of high-redshift cooling in 
the simulations does not have a dramatic effect on the cor- 
respondence between X-ray and total mass. We show there 
is a higher probability of finding an X-ray counterpart at 
high redshift, which can be attributed to a shorter time- 
scale on which ram pressure stripping could occur. By divid- 
ing our sample based on dynamical state we find subhaloes 
with Msub > 3 X lO^^/i~^M0 only lack an X-ray counterpart 
when the cluster is highly disturbed (in agreement with re- 
cent observations), however, relaxed clusters exhibit many 
deviations from the basic picture that light traces mass at 
lower subhalo masses. Joint weak lensing and X-ray analyses 
of relaxed systems are therefore also valuable and will yield 
much information about the physics of the ICM. 

A review of our main results in the presence of observa- 
tional noise and degraded resolution reveals many of these 
interesting mismatch scenarios are not currently observable, 
yet predicted improvements in lensing mass map resolution 
suggest these will be revealed in the coming decade, un- 
veiling frequent deviations from the simple assumption that 
light traces mass. 
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8 CONCLUSIONS 

In this paper, we have used resimulations of three non- 
radiative galaxy clusters in order to investigate the discrep- 
ancies between substructure in the hot gas and DM compo- 
nents, evident from recent comparisons of X-ray and weak 
lensing observations. We developed a simple technique to 
detect 2D substructures in simulated surface mass density 
and X-ray surface brightness maps of the clusters, without 
any reliance on circular symmetry or dynamical state. The 
resulting catalogues of 2D mass and 2D X-ray substructures 
were matched and we investigated how the success of this 
matching procedure varied with redshift, dynamical state 
and choice of gas physics employed. By utilising informa- 
tion about the underlying 3D subhalo distribution (obtained 
with SUBFINd) we have assigned subhaloes to the 2D mass 
substructures, allowing us to characterise the efficiency of 
our 2D substructure detection technique and reveal the ef- 
fect of subhalo mass on the 2D mass to 2D X-ray substruc- 
ture matching success. 

Our main results can be summarised as follows: 

• Having undertaken a thorough assessment of the prop- 
erties of the 2D substructure catalogues resulting from our 
novel detection procedure, we have ensured that any selec- 
tion effects or biases the technique may have introduced are 
understood. By attempting to match all 2D substructures 
detected in the surface mass density map with the 3D sub- 
haloes (identified with SUBFINd), we have concluded that our 
2D substructure catalogue is 90 per cent complete per mass 
bin (98 per cent overall) down to a 3D subhalo DM mass 
of ~ 10^^ ft- ^ Mo and 100 per cent complete down to a DM 
mass of 10"/i-^Mq. We are confident therefore that, in the 
3D subhalo mass range currently probed by weak lensing, 
the 2D substructure catalogues provide an accurate repre- 
sentation of the true 3D picture. We also establish that the 
2D mass substructure catalogue is pure and complete for 
^sub.TM > 10~^ Mpc, i.e. all 2D mass substructures 
with areas above this limit are successfully matched to a 3D 
subhalo and are, therefore, genuine. This purity threshold 
should allow the same detection procedure to be reliably 
applied to other simulated surface mass density maps in fu- 
ture, without the need for 3D subhalo data with which to 
compare. 

• We present a correlation between Asub.TM, the area of 
a 2D mass substructure, and Msub, the DM mass of the 
3D subhalo to which it is matched. The correlation is still 
apparent upon the introduction of basic observational noise, 
suggesting it could provide a quick estimate of the mass of 
a subhalo responsible for a peak in a weak lensing mass 
reconstruction, after accurate calibration. A measurement 
of the intrinsic scatter suggests such an estimate would be 
out by a factor of ~ 2. 

• The results of the matching between 2D mass substruc- 
tures and 2D X-ray substructures are surprising. We do not 
find X-ray counterparts for 23 — 33 per cent (depending on 
choice of X-ray catalogue) of all 2D mass substructures in 
the pure catalogue. Below Mgub ~ lO^"^/i~^M0 the match- 
ing success per mass bin begins to decrease significantly with 
decreasing subhalo mass. For the Icr X-ray catalogue, a few 
per cent of cluster cores, 5 per cent of group-size 2D mass 
substructures and 35 per cent of galaxy-size 2D mass sub- 
structures are not associated with a 2D X-ray substructure. 



The reasons for a matching failure are: 1) displacement of 
hot gas, where the X-ray substructure is intact yet spatially 
distinct from the DM, 2) depletion of hot gas, where so much 
gas has been stripped that detection of the 2D X-ray sub- 
structure fails or 3) complete disruption of the hot gas, where 
all hot gas appears to have been removed such that no 2D 
X-ray substructure is evident, even on visual inspection. We 
have conducted a detailed follow-up of examples of these 
scenarios with a set of case studies. 

• The dynamical state of the clusters (characterised by 
measuring the centroid shift variance in the X-ray surface 
brightness maps), is found to play a role in determining the 
fraction of 2D mass substructures without X-ray counter- 
parts. Substructures with Mgub > 3 x lO^'^/i~^M0 without 
an X-ray counterpart are restricted to the disturbed sample, 
suggesting major merger events are the cause. Substructures 
below this mass are less likely to have X-ray counterparts in 
relaxed systems, suggesting ram pressure stripping plays an 
important role on this scale; by definition, a long time has 
elapsed since the last merger, so these substructures have 
had most time to experience its effects. Similarly, the low 
redshift sample (0 ^ 2: ^ 0.2) contains more 2D mass sub- 
structures, in this mass range, that are unmatched to X-ray 
substructures than the high redshift sample (0.5 ^ z < 1). 

• The inclusion of high-redshift (until z ~ 5) cooling has 
only a mild impact on our results. It has little effect on the 
matching success between 2D mass substructures and 2D X- 
ray substructures for Msub > 3 x lO^'^/i~^M0, but reduces 
it, compared to the non-radiative case, below this due to the 
reduction of hot gas in these objects. 

• We have demonstrated that our simple 2D detection 
technique is still successful when noise which approximates 
that in real observations is added to the maps and the map 
resolution is degraded. As could be expected, the subhalo 
mass at which high completeness is achieved for the mass 
substructure catalogues is around an order of magnitude 
higher than in the fiducial set of maps. We have shown 
that this increase means many of the interesting mismatches 
which occur at lower mass scales cannot currently be ob- 
served. If the resolution of lensing mass maps can be im- 
proved by a factor of 10, to ~ 10/i~^kpc, we predict that 
many more discrepancies between the hot gas and dark mat- 
ter components of clusters will be observed and that these 
will not be restricted to rare, extreme merger events such 
as the Bullet cluster. Such an improvement, while dramatic, 
has been predicted for the coming decade and authors are 
already developing new observational analysis techniques to 
allow this, such that com parisons in the spirit of the present 
work can be undertaken l|Codl2009l ). These future observa- 
tions will provide a wealth of information about the physics 
of the ICM, the dynamical state of galaxy clusters and will 
allow us to probe the properties of the DM substructure 
directly. 



In future work, we will assess the impact of introducing 
heating processes, e.g. from galactic winds and the effects of 
active galactic nuclei into the resimulations, as well as the 
effect of including more realistic noise in the maps. 
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